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ABSTRACT 


This  is  Volume  2  of  the  final  report  summarizing  two  years  of  work  on  ground-based 
foliage-penetration  radar.  The  design  and  implementation  of  a  breadboard  radar  sys¬ 
tem  and  of  a  theoretical  and  experimental  investigation  of  target  and  propagation 
phenomena  specific  to  the  foliage -penetration  problem  are  described. 

The  radar  is  a  coherent  UHF  system  which  uses  a  large  (80-foot -diameter),  semi¬ 
circular  electronically  steered  array  antenna.  All  signal  processing  (including  pulse 
compression),  MTI  processing,  and  automatic  detection  use  digital  circuitry.  The 
digital  signal-processing  chain  ends  in  a  rather  sophisticated  thresholding  algorithm 
which  provides  constant -false -alarm-rate  detection  for  real-time  area  search  oper¬ 
ation  of  the  radar.  In  addition  to  this  area  search  mode,  the  radar  can  be  operated 
as  a  measurements  system  which  records  radar  measurements  digitally  on  magnetic 
tape  for  more  detailed  non-real-time  processing.  Finally,  the  radar  can  be  operated 
to  produce  real-time  periodograms  of  returns  from  a  single  range  azimuth  resolution 
cell.  This  mode  is  useful  for  monitoring  the  radar's  performance  as  well  as  for 
making  preliminary  investigations  of  various  target  and  clutter  returns.  Operating 
results  and  capabilities  are  described  briefly. 

Volume  1  comprises  three  main  sections:  a  short  introduction;  an  overall  description 
of  the  radar;  and  detailed  descriptions  of  the  hardware  and  software  subsystems  in 
the  radar. 

Volume  2  describes  the  use  of  the  radar  as  a  measurements  system  for  studying  tar¬ 
gets  and  the  clutter  environment,  in  addition  to  its  use  for  demonstrating  the  feasi¬ 
bility  of  radar  detection  of  walking  men  in  foliage.  Investigations  of  target  return 
spectra,  clutter  return  spectra,  and  clutter  return  amplitude  statistics  are  reported. 
Theoretical  models  relating  clutter  spectra  to  wind  turbulence  and  tree  resonances 
are  presented  and  compared  with  experimental  results.  Similarly,  a  previous  theo¬ 
retical  prediction  of  target  return  spectra  is  compared  with  experimental  results. 
Clutter  return  amplitude  statistics  are  measured  and  found  to  be  Gaussian  under  cer¬ 
tain  conditions.  Finally,  Vol.  2  contains  some  comments  on  the  relationships  between 
the  subclutter  visibility  and  various  parameters  of  digital  MTI  systems. 
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AN  EXPERIMENTAL  UHF  GROUND  SURVEILLANCE  RADAR 


ANALYSIS  AND  MEASUREMENTS 


I.  INTRODUCTION  R.  D.  Yates 

In  this  Volume,  the  spectra  of  clutter  and  target  returns  from  foliated  areas  are  presented, 
and  a  study  of  the  clutter  statistics  is  given.  In  addition,  an  expression  for  the  subclutter  visi¬ 
bility  in  digital  MT1  radars  is  presented. 

The  Long  Range  Demonstration  Radar  (LRDR)  program  was  terminated  on  30  June  1972. 

The  radar  as  a  detection  system  was  debugged  and  operated  satisfactor  ily  by  early  April  The 
measurements  system  part  of  the  radar  became  operational  by  the  middle  of  May.  Hence,  there 
was  a  limited  time  available  for  science  experiments. 

In  Sec.  II,  a  theoretical  model  is  presented  which  can  predict  the  clutter  spectrum  from 
foliated  areas  as  a  function  of  several  parameters.  Experiments  were  conducted  which  verified 
the  validity  of  the  model.  A  series  of  estimated  clutter  spectra  from  various  areas  under  dif¬ 
ferent  wind  conditions  is  presented. 

A  theoretical  model  which  can  predict  the  spectrum  of  a  target  moving  at  a  constant  velocity 
in  foliage  is  presented  in  Sec.  111.  The  estimated  spectra  of  a  transponder  used  as  a  point  target, 
and  those  of  a  single  man  are  presented.  It  appears  that  the  model  is  valid,  but  additional  tests 
arc  outlined  to  further  verify  the  model. 

In  Sec.  IV,  the  clutter  statistics  are  considered.  Based  on  theoretical  considerations  given, 
it  appears  that  the  temporal  statistics  may  be  modeled  as  a  Gaussian  random  process  over  the 
time  the  clutter  is  statistically  stationary.  An  experiment  is  presented  to  verify  this  model. 
Additional  tests  are  outlined  to  study  further  the  spatial  and  temporal  clutter  statistics. 

Finally,  Sec.  V  presents  an  expression  for  the  subcluttcr  visibility  in  digital  MTI  radars. 
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II.  CLUTTER  SPECTRUM 


A.  Theory  R.  D.  Yates 

M.  Labitt 

An  extensive  study  has  been  made  of  theoretical  and  experimental  models  for  the  clutter 

1  2 

power  spectral  density  from  wind-blown  foliage.  '  In  this  section,  the  results  will  be  recast  in 

radar  terminology  and  summarized. 

1 

Rosenbaum  and  Bowles  model  the  forest  by  a  flat  dielectric  slab  over  a  flat  conducting 
ground.  The  mean  or  background  slab  has  a  refractive  index  rj,  an  attenuation  factor  a,  and  a 
height  hj..  The  permittivity  of  the  slab  has  random  fluctuations  characterized  by  the  zero  mean, 
complex,  space-time  random  process  e(r,  t)  where  r  =  (x,  y,  z)  represents  the  transverse  direc¬ 
tion  {x,  y)  and  height  z,  and  t  is  the  time.  The  motion  of  the  scatterers  is  characterized  by  a 
velocity  field  V(  r,  t)  which  is  assumed  to  be  independent  of  the  permittivity  process  e(  r,  t).  The 
antenna,  which  is  horizontally  polarized,  is  assumed  to  be  high  enough  above  the  forest  so  that 
the  lateral  wave  propagation  mode  can  be  ignored. 

Although  general  expressions  for  the  clutter  power  spectral  density  are  derived,  the  most 
practical  situation  is  when  the  forest  is  assumed  to  be  statistically  homogeneous  and  isotropic 
in  the  transverse  directions  and  statistically  stationary  in  time. 

Under  the  above  assumptions,  and  using  a  range-gated  pulse-doppler  radar  such  as  the 
LRDR,  the  complex  envelope  of  the  clutter  signal  can  be  expressed  as 

X(t)  =  I(t)  +  jQ(t) 

_  G  1  2  •  p1/2  •  C(t)  (11-1) 

r  c 


in  which  I ( t )  and  Q(t)  are  the  in-phase  and  quadrature  components,  respectively,  and 
=  receiver  gain 

=  received  clutter  power  per  pulse  from  compressed  range  cell 

C(t)  =  unit  power,  complex,  clutter  random  process. 

The  received  clutter  power  may  be  written  as 


P  =  P  *  G  •  A  •  a  •  F 
CO  c  c 


(II-2) 


where 

Pq  =  transmitted  peak  power 

G  =  one-way  antenna  gain 

A  =  antenna  aperture  =  GA  /47r 

2 

a  =  clutter  cross  section  in  m 
c 

Fc  =  two-way  clutter  propagation  factor. 

Since  the  forest  is  assumed  to  be  statistically  homogeneous  and  isotropic  in  the  transverse 
directions,  the  clutter  cross  section  cr^  can  be  expressed  in  terms  of  an  area  cross  section  cr^  as 

a  =  a  •  R  *  AR  •  A9  (11-3) 

c  o 


2 


where 


2  2 

a  =  clutter  area  cross  section  in  m  /m 
o 


R  =  distance  to  range  cell 
AR  =  effective  range  cell  length  =  cT^/2 

Tp  =  effective  duration  of  the  compressed  pulse 
A0  =  two-way  antenna  beamwidth. 

On  letting  P(t)  denote  the  compressed  pulse,  its  effective  duration  is  defined  as 
/ 1  p(t)  | 2  dt 

t  =  — ; — rr  •  (II-4> 

p  max  |  P(t)  | 

Let  the  antenna  have  a  one-way  gain  pattern  P(0,  p)  where  G  and  </?  are  the  azimuth  and 
elevation  angles,  respectively.  Since  the  pattern  is  chosen  so  that  the  elevation  beamwidth  is 
large  and  ground  targets  are  of  prime  interest,  the  one-way  antenna  gain  G  and  two-way  beam- 
width  A0  are  as  follows: 


G  =  max  P(Q,  0)  =  P(0,  0) 
0 

JL  |P(O,0)|2dO 

A0  =  - - - p -  . 

|P(o,  o>r 


(H-5) 

<11— 6) 


Due  to  the  nature  of  the  clutter  model  used,  it  is  only  possible  to  express  the  product  cr^F^ 
in  terms  of  the  forest  parameters,  antenna  height,  and  distance  to  the  range  cell  as  follows: 


in  which 


a  F 
o  c 


47T 


2/2  .  2 
eo(7)  -1} 


•  R(0) 


h&  =  antenna  height 
h^.  =  slab  height 

A  =  free-space  wavelength 
R  =  range  to  cell 

=  free-space  permittivity 

T}  =  slab  refractive  index. 

The  function  R(r)  is  a  correlation  function  which  can  be  expressed  as 


(II-7) 


R(t)  =  f  f  \  f  <?( y-;  zi-  z2)  g(zi>  S’Hz2) 

X  m(z2,  t)  exp  [j  -yL  Jr]2  -  1  (z1  -  z 2)]  dz1  dz~ 


(II-8) 


3 


The  function  g(z)  represents  the  effects  of  ground  lobing  and  attenuation  within  the  forest  and 
is  expressed  as 

2 


in  which 


g(z)  =  exp  [-a  (l  -  ^)]  jl  -  exp  [-  ^  (a  +  jb>]  j 


2  a  77  hf 


47Thr 


(II”  9 ) 


(11-10) 


(II-H) 


The  function  m(z,r)  represents  the  effects  of  the  motion  of  the  scatterers  and  is  expressed 


m(z,  r)  =  E 


exp 


VT(z,t)  dt 


(H-12) 


where  V,p(z,  t)  is  the  component  of  the  velocity  field  V(x,  y,  z,  t)  in  the  transverse  direction  (x,  y), 
and  E{  }  is  the  statistical  expectation  operator.  If  a  carrier  at  (c/x)  Hz  is  frequency  modulated 
by  the  random  process  (47r/x)  .  V,p(z,  t),  then  m(z,  r)  is  the  correlation  function  of  the  modulated 
signal.  It  has  been  assumed  that  V(r,  t)  is  statistically  stationary  with  respect  to  t. 

The  function  <p{k;  z^,  z^)  represents  the  spatial  spectral  density  of  the  random  permittivity 
function  e(r,  t)  with  respect  to  the  transverse  directions  and  is  expressed  as  follows: 


<P( 


k;  Zy  z2)  =  ^  C(p;z1,z2)  exp  [- jp  •  k]  d2p 


<p(  |k|;z1,  z2)  due  to  statistical  homogeneity  and  isotropy 
in  the  transverse  directions. 


(11-13) 

The  function  C(p;  z^)  is  the  correlation  function  of  the  random  permittivity  and  is  defined  as 
C(p\zvz2)  -  E{t  (  R,  zy  t)  e*(  R  -  p,  zz>  t)} 


=  C ( |  p  j ;  Zy  z^)  due  to  statistical  homogeneity  and  isotropy 

in  the  transverse  directions,  and  stationarity 
in  time. 

The  clutter  process  C(r)  in  Eq.  (II-l)  has  a  correlation  function  RJt)  where 


(11-14) 


Rc(t)  =  E{C (t )  C*(t  -r)} 

.  R(t) 

"  R(0) 


(11-15) 


and  where  R(r)  is  given  by  Eq.  (II-8).  The  power  spectral  density  of  C(r)  is  denoted  as  Sc (f )  where 


Sc(f)  =  (  Rc(t)  e_;>27rfTdT  . 

Due  to  the  way  the  process  C(t)  has  been  defined,  it  possesses  unity  power. 


(11-16) 
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Although  it  is  not  possible  to  separate  the  product  three  propagation  terms  contained 

in  Eq.  ( II — 8 )  are  worth  noting.  If  the  random  fluctuations  of  the  permittivity  are  ignored,  the  two- 
way  propagation  factor  corresponding  to  a  point  at  height  h  within  the  slab  is  given  as 

,  I6r)2(h  -hf)4  , 

F(h)  =  - 2T  •  “1 - S - ~4  ■  lg(h)l  •  (11-17) 

(4ttR  )  (V  ~  1)  •  R 

The  first  factor  in  Eq.  (11-17)  is  the  free-space  loss,  the  second  is  the  loss  due  to  the  trans¬ 
mission  through  the  slab-free-space  interface,  and  the  third  is  the  loss  due  to  ground  lobing  and 
attenuation  within  the  slab. 

In  order  to  proceed  further  in  evaluating  the  clutter  spectrum,  the  statistics  of  the  velocity 
field  V(  r,t)  and  the  correlation  function  of  the  random  permittivity  C(p;  z^,  z^)  must  be  known. 
However,  there  are  two  limiting  cases  which  lead  to  simple  expressions  for  the  correlation  func¬ 
tion  m(z,  r) . 

Since  the  scatterers  are  bound,  the  velocity  process  V,p(z,t)  is  integrable  and  a  displacement 
process  D^(z,t)  can  be  defined  as 

Dx(z,t)=C  V  (z,  s)  ds  +  D  (z,  0)  .  (11-18) 

The  correlation  function  m(z,T)  can  then  be  written  as 

m(z,  t)  =  E^exp  { j  ^  [DT(z,  t)  -  L>T(z,  t  -  r) ]})  .  (11-19) 

A  phase  modulation  index  p(z)  can  then  be  defined  as 

M-(z)  =  ^  <j(z)  (11-20) 


ct2(z)  =  Var{DT(z,  t)} 


(11-21 ) 


The  first  case  of  interest  is  when  p(z)  «  1.  Then,  Eq.  (11-19)  can  be  approximated  as 

m(z,  t)  =  1  -  p2(z)  +  p2(z)  *  pd(z,t)  (11-22) 

in  which  p^(z,  t)  is  the  normalized  covariance  function  of  the  displacement  process  and  is  defined 
as 


Pd(z,t)  =  E 


DT(z,t)  -  DT(z,t) 
”  u(z) 


Dt(z,  t  -  t)  -  Drp ( z ,  t  -  t) 
ct(z) 


(11-23) 


In  order  to  derive  Eq.  (11-22),  the  statistics  of  Drp(z,t)  need  not  be  known.  Only  the  variance 

2  1 
a  (z)  and  normalized  covariance  function  p^(z,  r)  are  needed. 

The  second  case  of  interest  is  when  p(z)  »  1,  and  corresponds  to  the  case  of  large  index 

phase  modulation.  From  Ref. 3,  Eq.  (11-19)  can  be  approximated  as 


m(z,  r)  =  ^  (z, 

v 


(11-24) 


where  ^(z,  x)  is  the  characteristic  function  of  the  random  process  V^fz,  t)  and  is  defined  as 


5 


4'v(z,x)  =  E  {exp  [ jx  •  VT(z,t)]} 
=  ^  pv(z,  V)  eJxVdV 


(11-25) 


where  py(z,  V)  is  the  first-order  probability  density  function  of  the  process  V^,(z,  t). 

There  are  also  simplifications  of  Eq.  ( II— 8 )  for  the  correlation  function  R(r)  which  result 
from  assumptions  about  the  permittivity  process  c(r,  t).  The  first  case  is  when  it  is  assumed 
to  be  homogeneous  with  height.  Then,  the  statistics  of  the  velocity  field  do  not  depend  on  height, 
and  the  spatial  spectrum  function  ip( k;  z^,  z^)  is  a  function  of  z^  —  z^.  Thus,  Eq.  ( II  — 8)  reduces  to 


R(t)  =  R(0)  m(r) 

where 

R(0)  =  \  f\  f  z-i  -  z2)  g(zi)  g*(z2) 

x  exp|i  X"  -  1  <zi  _  z2*l  dzi  dz2 


(11-26) 


(11-27) 


m (t)  =  E  (exp  {j  y-  [D^,(t)  -  IX-ptt  -  t)]}) 


(11-28) 


It  follows  from  the  above  that  the  shape  of  the  clutter  spectrum  depends  only  on  the  statistics 
of  the  displacement  process  D^{t). 

The  second,  and  more  important,  simplification  of  Eq.  ( II  — 8)  is  when  the  spatial  spectrum 
function  can  be  written  as  follows: 

(p(k;  zy  z 2 )  =  V,2('X'  zi'  6(zi  -  z2*  (11-29) 


where  6(z)  is  the  delta  function.  Then,  Eq.  ( II— 8 )  can  be  expressed  as 

R(t)  =  R(0)  •  (  *  W(z)  m(z,  t)  dz 
^o 


(11-30) 


in  which 


W(z)  = 


</,2(T'’  z)  lg<z>!2 

fn  f  <P2(  Z)  |g(z)  |2dz 


(11-31) 


and 


R(0)  =  f  f  c pz(  z)  |  g(z)  |  2  dz 

r\ 


(11-32) 


Equation  (11-30)  would  result  if  one  assumed  that  the  scattering  was  uncorrclated  as  a  func¬ 
tion  of  height  and  that  the  effective  cross  section  was  proportional  to  W(z). 

2 

Based  on  heuristic  considerations,  Labitt  was  able  to  choose  a  weighting  function  W(z)  and 
a  correlation  function  m(z,r)  which  could  model  the  experimentally  measured  clutter  spectrum 
fairly  well.  This  model  is  described  below.  However,  the  function  R(0)  given  by  Eq.  (11-32), 
which  is  proportional  to  the  clutter  area  cross  section,  could  not  be  determined. 


6 


It  is  assumed  that  within  a  range-azimuth  cell,  all  the  trees  are  similar  and  that  their  mo¬ 
tion  may  be  approximated  by  a  linear  dynamic  system.  The  force  acting  on  the  trees  depends 
on  the  velocity  of  the  turbulent  wind  within  the  forest.  The  tree  displacement  then  phase  mod¬ 
ulates  the  transmitted  carrier  frequency. 

Let  the  average  wind  velocity  above  the  trees  be  denoted  as  v^  and  the  turbulent  component 
within  the  forest  as  v(t)  which  is  modeled  as  a  zero-mean,  statistically  stationary  random  proc¬ 
ess.  The  power  spectral  density  of  v(t)  may  be  approximated  as^ 


sv(f)  - 


_ k _ 

[i  +  (f/fe)2]5/6 


(11-33) 


where  k  is  a  constant,  and  f  is  a  cutoff  frequency.  The  power  of  the  wind  turbulence  depends 
in  some  manner  on  the  average  wind  velocity,  so  that  k  is  a  function  of  v  .  Typical  values  for 
the  cutoff  frequency  f<  appear  to  be  0.01  Hz. 

The  force  acting  on  the  trees  will  be  a  nonlinear  function  of  the  wind  turbulence  v(t).  Since 
the  trees  generally  have  a  highly  damped  motion,  it  would  seem  reasonable  to  assume  that  the 
foree  F(t)  eould  be  expressed  as 

F(t)  =  (v(t)]1"  v(t)  >  0 

=  —  l— v(t)J*'  v(t)  <  0  (11-34) 

where  v  is  a  small  positive  number.  It  is  also  reasonable  to  assume  that  v(t)  is  a  Guassian 
random  process.  Then  the  power  spectral  density  of  F(t)  can  be  expressed  as 


Sp(f)  =  £  AnS^n)(f)  (11-35) 

n  =  1 

where  S^n\f)  is  the  n^-fold  convolution  of  Sv(f),  and  are  decreasing  positive  constants.  Since 
S^(f)  is  impulsive,  except  for  the  exact  shape  about  zero  frequency,  Sp(f)  ean  be  approximated  as 


sF(f)  = 


[i  +  <f/f1)2]5/6 


(11-36) 


where  f^  is  some  convenient  cutoff  frequency,  and  a(v ■  )  determines  the  dependence  of  the  power 
in  the  force  F(t)  as  a  function  of  the  average  wind  velocity. 

The  tree  motion  is  modeled  by  a  linear  dynamic  system.  In  terms  of  an  electrical  analog, 
the  force  is  a  voltage  source  driving  a  series  RLC  cireuit  where  the  tree  displacement  is  pro¬ 
portional  to  the  charge  on  the  capacitor.  The  displacement  process  z,  t)  is  modeled  as 

Dt(z,  t)  =  <r(z)  *  d(t)  (11-37) 

where  a(z)  is  defined  by  Eq.  (11-21),  and  d(t)  is  a  normalized  displacement  process  with  unity 

power.  It  has  been  assumed  that  the  correlation  function  of  the  displacement  process  depends 

on  the  height  z  only  through  a  gain  faetor.  Sinee  trees  bend  like  a  beam,  the  height  dependence 

2 

of  the  standard  deviation  a(z)  was  taken  to  be  proportional  to  z  . 

Under  the  above  assumptions,  the  variance  and  normalized  covariance  functions  of  the  dis¬ 
placement  process  DT(z,t),  whieh  are  defined  by  Eqs.  (11-21)  and  (11-23),  respectively,  can  be 
written  as 
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(11-38) 


2  2  4 

ft  (z)  =  V  (VQ)  •  z 


PD(Z,  T)  =  pd(T) 


4 


Sd(f)  e^^df 


(11-39) 


in  which  /3(vq)  is  a  function  which  depends  on  the  average  wind  velocity  and  type  of  trees  and 
S^(f)  is  the  power  spectral  density  of  the  displacement  process  d(t)  and  is  expressed  as 


Sd(f) 


G  (f) 

V 


GT(f) 


where 


/Gv(f)  •  Grp(f)  df 
Gv(f)  =  [1  +  (f/fj)2]'5/6 


GT(f)  =  |1  + 


4,-1 


(11-40) 


(11-41) 


(11-42) 


In  Eq.  (11-41),  f^  is  the  cutoff  frequency  of  the  wind  turbulence;  in  Eq.  (11-42),  f  is  the  res¬ 
onant  frequency  of  the  trees  and  Qq  is  the  "QM  (or  quality)  factor.  Although  the  limiting  form  of 
the  correlation  function  m(z,r),  given  by  Eq.  (11-22),  is  now  known,  for  general  modulation  in¬ 
dices  the  statistics  of  d(t)  must  be  known.  It  was  decided  to  model  d(t)  as  a  Gaussian  random 
process  so  that  m(z,  r)  can  then  be  written  as 


m(z,  r)  =  exp  {—  n  (z)  [1  -  p  ,(t)]} 


where  p(z)  is  the  modulation  index,  and  from  Eqs.  (11-20)  and  (11-38)  is 
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/  \  47T 

p(z)  =  x 


0(vo) 


(11-43) 


(11-44) 


The  normalized  covariance  function  p^(r)  of  the  displacement  process  is  given  by  Eq.  (11-39). 

The  weighting  function  used  in  Eq.  (11-30)  must  now  be  determined  in  order  to  obtain  an  ex¬ 
pression  for  the  clutter  power  spectral  density  shape.  Again,  heuristic  considerations  were 
used  to  obtain  a  reasonable  answer. 

It  was  assumed  that  the  density  of  scatterers  as  a  function  of  height  is  given  as  y(z),  and 
that  the  field  intensity  I(z)  at  height  z  is  the  solution  of  the  equation: 

dl(z )  =  I (z )  y(z)  dz  (11-45) 


Hence,  I(z)  can  be  expressed  as 


I (z )  =  1(0)  ■  exp 


[£  Y(h)  ’  dh] 


It  follows  that  the  weighting  function  W(z)  can  then  be  expressed  as 

W(z)  =  I(z)  '  ?(z> 

'  ’  I(hf)  -  1(0)  ' 


(11-46) 


(11-47) 
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Several  choices  for  y(z)  were  tried  and  the  resulting  spectral  shapes  were  compared  with 

-  z2/2 

the  clutter  spectra  estimated  from  real  data.  The  best  choice  of  y(z)  was  found  to  be  e 
Thus,  the  weighting  function  can  be  written  as 


exp{-[L  •  erfc(z /\I2)  +  4-  z2]} 

W(z)  =  — ^ - - -  (11-48) 

f  exp  {— [  L  erfc(z /*J2)  +  z^]}  dz 


in  which  the  parameter  L  represents  the  total  propagation  attenuation  within  the  forest  and  is 
given  as 


L  =  In 


I(hf) 

Wi 


(11-49) 


and  erfc(x)  is  the  complementary  error  function.  As  a  practical  matter,  h^.  was  taken  to  be  in¬ 
finite  in  Eqs.  (11-48)  and  (11-49). 

Using  Eqs.  (11-39)  through  (11-44),  (11-48)  and  (11-30),  the  normalized  clutter  power  spectral 
density  given  by  Eqs.  (11-16)  and  (11-15)  ean  be  computed.  It  will  be  a  function  of  the  eutoff  fre¬ 
quency  of  the  wind  f^,  the  resonant  frequency  of  the  trees  fQ,  the  "Q"  of  the  trees  Qq,  a  standard 
deviation  gain  factor  /3(vq)  of  the  displacement  process  which  depends  on  the  average  wind  veloc¬ 
ity  v^  and  the  type  of  trees,  the  wavelength  A,  and  a  propagation  loss  factor  L  of  the  forest. 

One  could  question  some  of  the  assumptions  that  were  made  in  order  to  arrive  at  an  expres¬ 
sion  for  the  clutter  spectral  shape.  In  the  next  section,  experimental  verification  of  this  model 
will  be  presented. 


B.  Experiments 

1.  Clutter-Model  Verification  M.  I^abitt 

R.  D.  Yates 

The  clutter  model  developed  in  See.  A  above  was  tested  against  data  gathered  from  the  Florida 
Field  Site  Measurements  program.  The  radar  measurement  system  and  experiments  performed 
are  described  in  Ref.  5. 

The  quadrature  clutter  components  from  five  range  cells  were  recorded  simultaneously  at 
435  and  1  305  MHz.  The  range-azimuth  area  was  flat  land  with  pine  trees  and  light  underbrush. 

The  wind  velocity  above  the  range-cell  location  was  also  recorded. 

The  clutter  spectrum  was  estimated  at  U1IF  and  L-band.  The  tree  resonant  frequency  f  , 
tree  quality  factor  Qq,  displacement  standard  deviation  at  unit  height  /3(vq),  and  two-way  attenua¬ 
tion  L  were  chosen  so  that  the  theoretical  elutter  spectrum  matched  the  estimated  clutter  spec¬ 
trum  at  UlIF.  The  frequency  was  then  changed  to  L-band  and  the  attenuation  L  chosen  so  that 
the  theoretical  and  estimated  clutter  spectra  matched  up  again. 

The  displacement  spectrum  S^f),  and  the  corresponding  clutter  spectrum  S^lf)  at  43  5  and 
1  305  MHz,  based  on  the  theoretical  model,  are  shown  in  Fig.  1 1—  1 .  The  comparison  of  the  the¬ 
oretical  and  estimated  spectra  is  shown  in  Figs.  II-2  and  II  —  3  for  435  and  1305  MHz,  respectively. 

The  running  mean  and  standard  deviation  of  the  wind  velocity  above  the  trees  are  shown  in 
Fig.  II-4(a-b)  over  the  time  segment  used  to  estimate  the  clutter  spectrum.  It  is  seen  that  the 
wind  is  essentially  constant  at  lOmph. 
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The  parameters  used  in  the  theoretical  model  arc; 


UHF  L-Band 

Carrier  frequency  (MHz)  435  1305 

Tree  resonant  frequency  fQ  (Hz)  0.37  5  0.37  5 

Tree  quality  factor  Qq  3  3 

Tree  standard  deviation  at  unit 
height  £(v0)  (m)  0.021  0.021 

Two-way  attenuation  L  (dB)  15  40 

Wind  cutoff  frequency  fc  (Hz)  0.01  0.01 


Other  clutter  areas  were  considered  and  similar  results  were  obtained.  Based  on  this,  it 
would  appear  that  the  clutter  model  presented  in  Sec.  A  can  be  used  to  predict  the  shape  of  the 
AC  clutter  spectrum  as  a  function  of  the  parameters  f  ,  f^,  Q^,  /3(v /,  X  and  L. 

It  was  not  determined  whether  the  exact  shape  of  the  clutter  spectrum  over  the  range  of 
±0.1  Hz  was  predicted  correctly.  The  wind  cutoff  frequency  f  did  not  affect  the  shape  of  the  clut¬ 
ter  spectrum  outside  this  region.  However,  it  greatly  affects  the  apparent  DC  power  present. 

The  shape  of  the  scatterer  density  *y(z)  about  z  =  0,  as  given  in  Eq.  (11-45),  will  affect  the  amount 
of  DC  power  present  also.  When  a  clutter  spectrum  was  estimated  with  the  DC  left  in,  the  differ¬ 
ence  between  the  theoretical  and  estimated  DC  clutter  power  was  at  most  10  dB. 

It  should  be  noted  that,  given  a  particular  range  azimuth  cell,  the  observed  DC  clutter  power 
will  be  a  random  variable.  For  a  homogeneous  forest,  the  DC  clutter  power  will  be  exponentially 
distributed  over  the  collection  of  range  cells.  This  makes  it  difficult  to  match  up  the  theoretical 
and  estimated  clutter  spectra  near  zero  frequency. 

2.  Tree  and  Wind  Experiments  M.  Eabitt 

R.  D.  Yates 

In  the  clutter  model  developed  in  Sec.  A,  the  phase  modulation  index  given  by  Eq.  (11-44)  is 

4  TT  2 

H(z)  =  —  •  /3(vq)  •  z  .  (11-50) 

The  function  £(v^),  where  vq  is  the  average  wind  velocity  above  the  trees,  can  be  interpreted  as 

the  standard  deviation  of  the  displacement  process  D(z,t)  at  unit  height.  As  vq  increases,  the 

wind  turbulence  within  the  trees  also  increases,  so  that  one  could  postulate  that  0(vq)  is  a  mon- 

otonically  increasing  function  of  v  .  When  £(vq)  is  known,  the  behavior  of  the  clutter  spectrum 

as  a  function  of  v  can  then  be  determined, 
o 

In  the  clutter  model,  a  displacement  process  D(z,  t)  was  defined  and,  based  on  heuristic  ar¬ 
guments,  an  expression  for  its  power  spectrum  was  given  by  Eq.  (11-40).  As  seen  from  Sec.  1 
above,  the  theoretical  and  estimated  clutter  spectra  could  be  made  to  agree.  However,  the  dis¬ 
placement  spectrum  could  only  be  measured  indirectly  for  low  modulation  indices. 

An  experiment  designed  to  check  the  validity  of  the  displacement  spectrum  and  to  estimate 
the  function  /3(vq)  is  depicted  in  Fig.  II - 5 .  A  thin  horizontal  wire  was  attached  to  a  tree  trunk 
3m  above  the  ground,  looped  over  a  stationary  pulley,  and  attached  to  a  small  weight.  The  pul¬ 
ley  shaft  was  connected  to  a  precision  servo  potentiometer.  A  DC  voltage  was  applied  to  the  po¬ 
tentiometer,  the  output  of  which  was  connected  to  an  A/D  converter.  A  hot  wire  anemometer 
was  used  to  measure  the  wind  speed  at  the  treetop  level  in  the  area  of  the  tree.  Its  output  was 
also  connected  to  an  A/D  converter.  The  digitized  data  were  recorded  on  magnetic  tape  in  the 
same  manner  as  for  the  clutter  and  target  data. 
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The  power  spectral  density  estimate  of  the  wind  velocity  is  given  in  Fig.  11-6.  As  expected, 
an  f  ^  roll-off  in  frequency  was  observed. 

The  power  spectral  density  estimate  of  the  tree  displacement  is  given  in  Fig.  11-7.  In  addi¬ 
tion,  the  theoretical  displacement  spectrum  given  by  Eq.  (11-40)  is  shown. 

In  estimating  the  displacement  spectrum,  the  DC  component  was  removed,  hence,  the  notch 
at  zero  frequency.  It  is  seen  that  the  estimated  displacement  spectrum  is  wider  about  ±0.1  Hz 
than  the  theoretical  spectrum.  There  are  also  two  spikes  present  at  ±0.9  Hz. 

With  the  anemometer  30  m  from  the  tree,  the  squared  coherency  function  between  the  wind 
velocity  and  tree  displacement  was  estimated,  as  shown  in  Fig.  II-8.  One  would  expeet  no  cor¬ 
relation  between  the  tree  displacement  and  wind  velocity.  It  is  seen  that  a  correlation  exists  at 
±0.91Iz  and  about  zero  frequency.  Based  on  these  results,  it  was  dedueed  that  the  spectral  broad¬ 
ening  about  0  Hz  in  the  displacement  spectrum  was  due  to  system  instabilities.  The  lines  at 
±0.9IIz  were  due  to  a  resonance  in  the  measurement  system,  possibly  in  the  wire  or  potentiometer. 

Based  on  the  above  results,  one  can  conclude  that  the  theoretical  expression  for  the  displace¬ 
ment  spectrum  is  accurate. 

The  next  problem  considered  was  the  behavior  of  the  function  £(v^).  Due  to  the  fact  that  the 
displacement  spectrum  has  most  of  its  power  near  zero  frequency,  a  measurement  of  its  standard 
deviation  would  be  dominated  by  low-frequency  components.  The  displacement  process  model 
assumed  that  the  displacement  spectrum  shape  was  independent  of  the  average  wind  speed.  It 
follows  that  the  variance  of  the  derivative  of  the  displacement  process  can  be  expressed  as 

=  /3Z(vo)  •  z4  •  |/>u(°)|  (11-51) 

where  p^(o)  is  the  second  derivative  of  the  normalized  displacement  covariance  function.  It  fol¬ 
lows  that  the  behavior  of  /3(vq)  can  also  be  determined  from 

The  tree  displacement  samples  and  corresponding  wind  velocity  samples  were  divided  into 
eonseeutive  bloeks  of  6  minutes  in  duration.  The  sample  mean  of  the  wind,  and  sample  var¬ 
iance  tx  ?  of  the  tree  velocity  were  computed  for  each  block.  The  plot  of  vs  vq  is  shown  in 
Fig.  II-9.  A  least-mean-squarc-error  fit  of  cr^  =  a  •  v^  indicated  that  b  =  1.47,  with  a  correla¬ 
tion  coefficient  of  0.937.  The  estimate  of  the  variance  of  the  tree  velocity  was  computed  by  tak¬ 
ing  discrete  Fourier  transforms  over  the  block,  multiplying  the  magnitude  squared  of  the  fre- 

2 

quency  coefficients  by  f  ,  and  summing  all  the  modified  coefficients. 

Based  on  these  results,  it  is  postulated  that  the  function  /3(vq)  has  the  following  form: 

0<vo)  =  K  ■  vo15  (11-52) 

where  K  depends  on  the  type  of  forest. 

It  should  be  remarked  that  Eq.  (11-52)  was  deduced  for  wind  speeds  between  0  and  20rnph. 

It  is  not  clear  if  this  expression  is  valid  well  beyond  20mph. 

3.  Clutter  Speetrum  Estimates 

The  measurement  mode  of  the  LRDU  was  used  to  record  the  quadrature  components  of  the 
elutter  return  from  five  eonseeutive  range  cells  from  three  foliage  areas  under  different  wind 
conditions.  The  wind  velocity  at  the  top  of  the  antenna  tower  was  also  recorded. 
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(a)  Spectral  Estimation  Scheme 


M.  Labitt 
R.  D.  Yates 


A  detailed  discussion  of  clutter  and  target  returns  and  the  estimation  of  their  power  spec¬ 
tral  density  is  given  in  Ref.  5.  A  summary  of  the  estimation  scheme  is  presented  here  for 
completeness. 

In  the  LRDR,  the  compressed,  range-gated,  quadrature  video  signal  samples  are  produced 
at  a  rate  of  (555.6  X  10  S  *  =  1799.8560Hz  per  channel.  The  TAPWRITE  program  reduces  the 
data  rate  to  112.4910  Hz  by  segmenting  the  data  into  consecutive  blocks,  each  block  containing 
16  points.  These  points  arc  added  together  to  form  one  new  data  point.  The  quadrature  compo¬ 
nents  from  five  consecutive  range  cells  are  then  recorded  on  digital  magnetic  tape. 

A  series  of  computer  programs,  written  to  analyze  the  data  from  the  Florida  Everglades 
measurement  program,  was  used  to  study  the  data  produced  by  the  LRDR  measurement  program. 
The  basic  procedure  may  be  divided  into  three  steps. 

The  first  step  is  to  filter  the  data  further,  depending  on  the  frequency  extent  of  the  signals 
being  studied.  The  data  per  quadrature  channel  were  divided  into  consecutive  groups  of  N&  points 
per  group.  The  points  in  each  group  were  then  added  to  produce  one  filtered  data  point.  The 
sampling  rate  of  the  filtered  data  is  then 


f 

s 


112.4910 

X- ~~ 

a 


IIz 


(11-53) 


The  second  step  consists  of  forming  a  set  of  frequency  coefficients  using  the  FFT  (fast 
Fourier  transform)  algorithm.  The  data  per  quadrature  channel  are  divided  into  B  consecutive 

blocks,  each  block  containing  N  points.  Let  ^  *  be  the  N  samples  of  the  in-phase 

th  n"°  (b) 

and  quadrature  components  of  the  bLn  block.  Let  the  complex  signal  sample  Xn  ]  be  defined  as 


X<b)  =  i(b)  (b) 

n  n  J  n 


(11-54) 


The  N  frequency  coefficients 


N-l 

k=0 


for  the  b**1 


block  are  computed  as 


(b)  _ 


N-l 

Z 

n=0 


X 


(b) 


W 


r  .27 r 

exp [-3  ^ 


nkl 


(11-55) 


in  which  { W^}  are  a  set  of  time  weights  used  to  obtain  the  appropriate  frequency  resolution.  In 
all  cases,  they  were  chosen  as 


W 

n 


.  2  7T  ,  ^  1  . 

sin  (n  +  y ) 


(11-56) 


In  Eq.  (11-55),  N  was  constrained  to  be  a  power  of  2.  The  spectral  estimates,  based  on  one  FFT, 
are  then  defined  as 


pr-iA^i2 


(11-57) 


The  third  step  is  to  smooth  the  estimates  given  by  Eq.  (11-57)  using  the  data  from  the  B  blocks. 

||  M-l 

P,  I  ,  were  computed  as 
k,k  =  0 
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(11-58) 


B 


Pk  “  B  -  Pk 
b  =  l 

(b)  . 


(b) 


The  samples  X n  in  Eq.  (11-55)  can  be  interpreted  as  a  continuous  signal  X(t)  sampled  every 


l/f  sec.  Hence, 
'  s 


■(b)  . 


=  X  |(n  +  bN)  — 


(11-59) 


On  assuming  that  the  signal  X(t)  is  a  stationary  random  process  with  a  power  spectral  den¬ 
sity  Sx(f),  and  that  the  sampling  is  high  enough  so  that  no  aliasing  takes  place,  the  estimates  of 
the  spectrum  S^(f)  are  then  taken  as 


N-k 


for  k  =  0,  1,  2,  . 

for  k  =  —  1,  —  2, 


N 

2 


N  -  2 


The  expected  value  of  can  be  shown  to  be 

#-■) 


pfs/2  _  /kf„ 

E{Pk}=  \_S 


S  (f)  R 
-(f  /2)  X 


df 


(11-60) 


(11-61) 


in  which  R(f)  is  a  frequency  resolution  function,  and  Sx(f)  is  the  aliased  version  of  S^( f) . 
The  resolution  function  R(f)  is  given  in  terms  of  the  time  weights  {Wn)  as 


R(f)  = 


N-1 

V  ,r  27rnf 

L  WnexP-J  — 

n=0  & 


(11-62) 


The  aliased  spectrum  Sx ( f )  is  expressed  in  terms  of  S^(f)  as 


OO 

S  (f)  =  V  S  (f-mf  ) 
x  x  s 

m=-°° 


(11-63) 


On  assuming  that  the  spectral  estimates  from  block-to-block  are  uncorrelated,  the  variance 
of  is  given  as 

Var  {P  (b)} 

Var{Pk}  =  - g-* -  .  (11-64) 


The  right-hand  side  of  Eq.  (11-64)  does  not  depend  on  b  due  to  the  assumption  that  X(t)  is  statis¬ 
tically  stationary.  If  X(t)  is  a  Gaussian  process,  then  Var{P^^}  =  E2{P^^},  so  that 


Var{Pk}  * 
E2{Pk}  K 


(11-65) 


The  spectral  estimate  of  a  random  process  based  on  one  FFT  is  seen  to  have  a  variance  governed 
by  Eq.  (11-64)  with  B  =  1.  Hence,  in  order  to  obtain  a  reliable  estimate,  B  should  be  large.  It 
was  found  that  adequate  estimates  of  the  clutter  spectrum  could  be  obtained  with  B  in  the  range 
of  50  to  100. 
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The  resolution  of  the  estimates  is  governed  by  the  choice  of  {\Vn}. 
gular  bandwidth  of  R(f)>  over  its  main  period  |f|^  fg/2,  is  defined  as 

rfs/2 

/.< fs/2)  R(f)  df 

k  ~  max  R(f) 

|fk(fsA) 


The  equivalent  rectan- 


(11-66) 


On  letting  the  weights  {Wn}  be  real,  and  normalizing  R(f)  so  it  has  unit  area  over  its  main  per¬ 
iod,  it  follows  that 


pf  /2 

J-(fs/2) 


R(f)  df  =  1 


(11-67) 


R(— f)  =  R(f)  (11-68) 

max  R(f)  =  R(0)  =  .  (11-69) 

|fk<fs/2)  r 


On  using  the  weights  given  by  Eq.  (11-56),  it  follows  that  the  equivalent  rectangular  bandwidth 
is  given  as 


f  = 
r 


(11-70) 


The  resolution  function  has  zeros  at  f  =  ±kfg/N  for  k  =  2,  3,  .  .  . ,  N/2,  and  the  maximum  of  the 
sidelobes  occurs  approximately  at 


f  =  ± 


(k  +  j)  f 
N 


s 


for  k  =  2,  3,  .  .  .  ,  (N  —  2)/2.  The  first  6  sidelobes  are  31.5,  41.5,  48.5,  54,  58,  and  62  dB  below 
the  main  peak.  The  value  of  R(f)  at  ±  fg/N  is  6  dB  below  the  main  peak.  The  equivalent  rectan¬ 
gular  bandwidth  f  can  be  used  as  a  measure  of  the  resolution  of  the  spectral  estimates. 

Given  that  R(f)  possesses  adequate  resolution,  it  follows  that  Eq.  (11-61)  can  be  approximated 
as 


/kf  \  /kf  \ 

E<pk>  =  pDc-  r\w)  +  sac(-w)  <n-71> 

where  it  has  been  assumed  that 

Sx(f)  =  PDC6(f)  +  SAC(f)  *  (II'72) 

Hence,  if  the  spectrum  of  X(t)  has  a  DC  component,  it  will  only  appear  in  the  estimates 

P  P  P 
*0'  1'  -1* 

In  practice,  the  statistics  of  the  process  X(t)  may  not  be  stationary  across  the  B  blocks  used 
to  compute  the  spectral  estimates.  In  the  case  of  clutter,  when  the  phase  modulation  index  is  less 
than  1,  it  follows  from  Sec.  II-A  that  the  AC  clutter  spectrum  shape  does  not  vary  as  a  function  of 
the  average  wind  velocity.  However,  the  AC  clutter  power  does.  By  estimating  the  AC  power  in 

each  block,  and  normalizing  the  data  to  unit  AC  power,  the  smoothing  from  block-to-block  should 

2  th 

be  more  efficient.  The  sample  mean  and  sample  variance  ,  for  the  bLn  block,  were  com¬ 
puted  as  follows: 
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(11-73) 


N-l 


Mb  = 


N 


■(b) 


n  =  0 


N-l 
1  V 


N-l  U 
n= 0 


|X(b)  -M,  |2 
1  n  b 1 


(11-74) 


The  normalized  data  from  the  bth  block  were  then  taken  as 


,M  X(b)-Mh 
!  (b)  _  n  b 


(11-75) 


The  spectral  estimates  based  on  the  normalized  samples  are  then  computed. 

It  was  also  of  interest  to  estimate  the  cross  correlation  between  different  signals.  Let  the 
two  stationary  random  processes  X(t)  and  Y(t)  have  power  spectra  S  (f)  and  S  (f)  and  a  cross 

^  y 

spectrum  Sx^(f).  The  function  of  interest  is  the  squared  coherency  function  defined  as 


C  (f)  =  c  m 
xy'  S  (f) 


Vf)| 


s  (f) 

y 


(11-76) 


The  cross  spectrum  S  (f)  is  defined  as  the  Fourier  transform  of  the  cross-correlation  function 

xy 

R  (t)  =  E{X(t)  Y*(t  —  t)}.  If  there  is  no  correlation  between  X(t)  and  Y(t),  then  C  (f)  =  0.  On 
xy  xy 

the  other  hand,  complete  correlation  leads  to  C  (f)  =  1.  It  has  been  assumed  that  X(t)  and  Y(t) 

xy 

contain  no  specular  components  so  that  no  impulses  are  present  in  the  spectra. 

The  squared  coherency  function  was  estimated  by  assuming  the  spectra  of  X(t)  and  Y(t)  as 
already  described,  and  the  cross  spectra  by  replacing  |A/b^  |  in  Eq.  (11-57)  by  A  ^  •  A  ^  , 

k  x,k  y,k 

where  *  denotes  the  complex  conjugate  operation,  A  /  the  frequency  coefficients  of  the  process 
(b)  X'K 

X(t),  and  A  ^  those  of  Y(t).  The  estimate  of  the  squared  coherency  function  was  then  taken  as 


C  (f,  ) 
xy  k 


S  (f.  ) 

1  xy  k  1 


Sx<fk> 


S  (f.  ) 
y  k' 


(II-77) 


On  assuming  adequate  resolution  and  the  absence  of  specular  components,  it  can  be  shown 


that 


E{ C  (f,  )}  =  C  (f,  )  +  4 
k  xy  k  ’  xy  k  B 


(11-78) 


Var{C  (f.  )}  =  4  [C  (f.  )  +  C2  (f.  )]  +  -k 
1  xy  k  J  B  1  xyv  k'  xy  k  J  ^2 


(H-79) 


where 


fk  = 


kf 

_ s 

N 


(11-80) 


The  estimation  schemes  described  were  extensively  tested  against  both  real  and  simulated 
data.  It  was  found  that  the  power  spectral  densities  and  coherency  functions  could  adequately 
be  estimated  by  the  methods  described  in  this  section. 
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(b)  Results 


R.  D.  Yates 


The  estimated  clutter  spectra  and  squared  coherency  functions  are  shown  in  Figs.  11-10 
through  11-29.  Three  areas  were  considered,  each  one  consisting  of  five  consecutive  range  cells. 
The  range  cells  were  1  5  m  in  length  and  the  azimuth  beamwidth  was  2.5°.  The  parameters  of 
range-cell  areas  studied  are  given  in  Table  II- 1. 


TABLE  11-1 

PARAMETERS  OF  RANGE-CELL  AREAS  STUDIED 

Area  1 

Area  2 

Area  2 

Area  3 

Area  3 

First  range-cell 
number 

313 

413 

413 

1120 

1120 

Beam  number 

26 

5 

5 

26 

26 

Range  to  first  cell 
(km) 

5.7 

7.2 

7.2 

17.8 

17.8 

Approximate  wind 
speed  at  antenna 
tower  (mph) 

5  to  10 

10  to  2  0 

10  to  20 

5  to  10 

10  to  20 

Tape  number 

28 

31 

31 

29 

22 

Records  analyzed 

1006  to  2005 

401  to  800 

11  to  400 

6  to  2400 

6  to  2005 

30-dB  bandwidth 
of  clutter  spectra 
(Hz) 

±1.43 

±2.81 

±3.84 

±1.48 

±2.25 

Figure  numbers 

11-10  to  -16 

11-17  to  -19 

11-20  to  -22 

11-23  to  -26 

11-27  to  -29 

The  estimates  of  the  clutter  power  spectrum  for  each  range  cell  are  presented  on  a  decibel 
vs  frequency,  and  a  decibel  vs  log-frequency  scale.  The  latter  case  is  useful  in  determining  the 
slope  of  the  clutter  spectrum  tails.  The  average  of  the  spectra  of  the  five  range  cells  is  also 
given.  The  squared  coherency  function  between  selected  range  cells  is  given  to  determine  if  any 
correlation  exists.  If  there  were  antenna  motion  for  example,  the  DC  clutter  in  each  range  cell 
would  be  modulated  in  an  identical  manner.  This  common  signal  would  then  appear  in  the  coher¬ 
ency  function. 

In  most  tests,  the  wind  speed  at  the  top  of  the  antenna  tower  was  recorded.  The  running 
mean  and  standard  deviation  of  the  wind  speed  for  the  time  duration  of  the  data  used  to  estimate 
the  spectra  are  presented.  Since  the  range  cells  were  at  least  5.7  km  from  the  radar,  the  wind 
speed  at  the  range-cell  location  could  be  different.  If  the  clutter  spectrum  bandwidth  is  to  be 
correlated  with  the  wind  velocity,  a  wind  sensor  should  be  placed  at  the  range-cell  location. 

Area  1  was  5.7  km  from  the  radar.  The  clutter  spectral  estimates,  with  the  DC  filtered  out 
before  the  estimation  was  made,  are  shown  in  Figs.  11-10  through  11-13;  Figs.  11-14  through  11-16 
show  the  estimated  spectra  of  the  same  data  with  the  DC  left  in. 

It  is  seen  that  the  spectra  have  the  shape  predicted  by  the  clutter  model  presented  in 
Sec.  II-A-1.  The  resonant  frequency  of  the  trees  is  seen  to  be  0.336Hz.  With  the  DC  removed, 
it  appears  that  the  shape  of  the  spectra  in  the  region  ±0.140  Hz  is  identical  to  the  turbulent  wind 
velocity  spectra.  The  tails  of  the  clutter  spectra  have  a  slope  slightly  less  than  —  5§  in  cells 
1  through  4.  In  cell  5,  the  slope  is  —  5f .  The  flat  section  of  the  spectra  past  ±3  Hz  is  due  to  the 
receiver  noise.  It  is  seen  that  the  DC  clutter  power  from  cell-to-cell  in  Figs.  II— 1 4(b)  through 
(f)  is  not  constant.  The  width  of  the  spectra  30  dB  below  the  tree  resonant  peaks  is  ±1.43  Hz. 
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System  tests  indicated  that  the  DC  offsets  in  the  quadrature  video  channels  were  negligible. 
Hence,  the  DC  powers  shown  are  entirely  due  to  the  clutter  signal.  When  comparing  the  DC-to- 
AC  power,  it  should  be  remembered  that  the  estimated  DC  power  at  0  Hz  is  equal  to  K  •  •  f 

where  K  is  a  system  constant,  is  the  DC  clutter  power,  and  f^  is  the  frequency  resolution. 

The  AC  spectrum  is  equal  to  K  •  S^^(f),  where  S^,(f)  is  the  AC  clutter  spectrum,  when  adequate 
resolution  exists.  Tests  were  made  to  ensure  that  adequate  resolution  was  used  in  all  the  spec¬ 
tral  estimates  presented. 

The  squared  coherency  function  between  several  range-cell  pairs  is  shown  in  Figs.  II- 12(a) 
through  (e)  with  the  DC  absent,  and  in  Figs.  II- 16(a)  through  (e)  with  DC  present.  If  no  correla¬ 
tion  exists  between  range  cells,  then  the  mean  and  standard  deviation  of  the  estimate  of  the  squared 
coherency  function  are  B  where  B  is  the  numher  of  data  blocks  used  to  smooth  the  estimates. 

It  is  seen  that  no  correlation  exists  between  nonadjacent  range  cells.  There  is  partial  correlation 
between  adjacent  cells,  as  seen  in  Figs.  II-12(a)  and  (e),  due  to  the  fact  that  the  basic  100-nsec 
pulse  is  passed  through  a  10- MHz  filter  at  the  receiver  which  causes  the  apparent  range-cell 
length  to  be  wider  than  15  m. 

Area  2  was  7.2  km  from  the  radar.  Spectral  estimates  for  two  different  wind  conditions  are 
shown  in  Figs.  11-17  through  II-  19  and  11-20  through  11-22,  respectively.  The  frequency  resolution 
used  is  identical  in  both  cases.  However,  due  to  the  extent  of  the  spectrum,  the  tree  resonant 
frequency  could  not  be  resolved. 

In  the  first  case,  the  30-dB  bandwidth  is  ±2.81  Hz,  while  for  the  second,  it  is  ±3.84  IIz.  These 
were  the  widest  clutter  spectra  seen  thus  far.  Comparing  corresponding  range  cells  in  both  cases, 
it  is  seen  that  as  the  AC  clutter  spectrum  spreads  out,  the  DC  power  is  reduced.  Again,  this  was 
predicted  bv  the  clutter  model. 

Area  3  was  18  km  from  the  radar.  Spectral  estimates  for  two  different  wind  conditions  are 
shown  in  Figs.  11-23  through  11-26  and  11-27  through  11-29,  respectively.  The  first  case  is  for  a 
low  wrind  velocity.  The  clutter  spectrum  has  the  shape  as  predicted  by  the  clutter  model.  The 
tree  resonant  frequency  in  this  case  is  0.  3  Hz.  The  slope  of  the  clutter  spectrum  tails  is  seen  to 
be  approximately  -  5f .  There  is  no  correlation  betw  een  nonadjacent  range  cells.  The  30-dB 
bandwidth  is  seen  to  be  ±1.48  Hz.  Within  the  experimental  error,  this  is  the  same  as  the  band¬ 
width  for  the  clutter  spectra  from  area  1. 

The  second  case  is  for  a  higher  wind  velocity.  The  30-dB  bandwidth  has  increased  to  ±2.25. 
Comparing  corresponding  range  cells  for  the  two  cases,  it  is  seen  that  as  the  AC  clutter  spectrum 
spreads  out,  the  DC  clutter  power  is  reduced.  Again,  there  is  no  correlation  between  nonadjacent 
range  cells. 
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Fig.  II- 1.  Theoretically  computed  displacement  Fig.  II-2.  Comparison  between  theoretical 
spectrum  S^(f),  and  clutter  spectra  S^(f ).  and  estimated  clutter  spectra  at  4  35  MHz. 


Fig.  II- 3.  Comparison  between  theoretical  and  estimated 
clutter  spectra  at  1305  MHz. 
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MILES  PER  HOUR 
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TIME  (min.) 

(b) 

Fig.  II- 4.  Wind  data  for  clutter  tape:  (a)  running  standard  deviation  of  wind; 
(b)  running  average  of  wind. 
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Fig.  II- 5.  Wind  and  tree  motion  instrumentation. 
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Fig.  II- 6 .  Wind  spectral  estimation.  A  hot  wire  anemometer 
was  used  as  sensor.  As  expected,  curve  approximates  an  f_5/3 
roll-off  observed  by  many  other  experimenters. 


Fig.  II- 7.  Tree  displacement  spectrum.  Spikes  at  ±0.9  Hz  are  spurious. 
Leveling  off  at  —  65dB  caused  by  system  noise.  Notch  at  0  Hz  because 
algorithm  removes  DC. 
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Fig.  II—  8.  Squared  coherency  spectral  estimation  between  wind  and  tree 
displacement.  Here,  wind  and  tree  sensors  were  more  than  30  m  apart. 
One  would  not  expect  any  correlation,  but  some  is  observed  near  0  and 
0.9  Hz;  these  are  believed  to  be  spurious. 


Fig.  II- 9 .  Estimated  standard  deviation  of  tree  velocity  vs  estimated 
wind  speed.  Superimposed  is  a  least-mean-square-error  power-law 

fit.  o-^  =  a  •  v  F  47. 

D  o 
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RELATIVE  POWER  ( D  B  >  RELATIVE  POWER  (DB>  RELATIVE  POWER  ( DB ) 


FREQUENCY  i HZ t 
SPECTRAL  ESTIMATION 


FREQUENCY  (H 7) 
SPECTRA!,  FSTIMATION 


Fig.  II- 1 0.  Clutter  from  area  1.  DFT  length  =  256,  resolution  =  0.0659  Hz,  sampling 
rate  =  11.25  Hz.  (a)  DFTs  averaged  270;  (b)  through  (f)  DFTs  averaged  =  54. 
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Fig.  11-11.  Clutter  from  area  4.  DFT  length  =  256,  resolution  =  0.0659  Hz,  sampling 
rate  =  11.25  Hz.  (a)  DFTs  averaged  =  270;  (b)  through  (f)  DFTs  averaged  =  54. 
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Fig.  II-12{a-e).  Clutter  from  area  1;  correlation  between  range  cells. 
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Fig.  11-13.  Clutter  from  area  1;  wind  speed  at  antenna  tower,  (a)  Running  standard 
deviation  of  wind;  (b)  running  average  of  wind. 


25 


RELATIVE  POWER  ( DB >  RELATIVE  POWER  (DB)  RELATIVE  POWER  (DB) 


FREQUENCY  (HZ) 
SPECTRAL  ESTIMATION 


FREQUENCY  (HZ) 
SPECTRAL  ESTIMATION 


Fig.  11-14.  Clutter  from  area  1.  DFT  length  =  256,  resolution  =  0.0659  Hz,  sampling 
rate  =  11.25  Hz.  (a)  DFTs  averaged  =  270;  (b)  through  (f)  DFTs  averaged  =  54. 
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Fig.  II- 1 5.  Clutter  from  area  1.  DFT  length  256,  resolution  =  0.0659  IIz,  sampling 
rate  -  11.25  Hz.  (a)  DFTs  averaged  -  270;  (b)  through  (f)  DFTs  averaged  =  54. 
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Fig.  II-l6(a-e).  Clutter  from  area  1;  correlation  between  range  cells. 
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Fig.  11-17.  Clutter  from  area  2.  DFT  length  256,  resolution  =  0.165  Hz,  sampling 
rate  28.12  Hz.  (a)  DFTs  averaged  -  270;  (b)  through  (f)  DFTs  averaged  -  54. 
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Fig.  11-18.  Clutter  from  area  2.  DFT  length  =  256,  resolution  =  0.165  Hz,  sampling 
rate  =  28.12  Hz.  (a)  DFTs  averaged  =  270;  (b)  through  (f)  DFTs  averaged  =  54. 
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Fig.  11-19.  Clutter  from  area  2;  wind  speed  at  antenna  tower,  (a)  Running  standard 
deviation  of  wind;  (b)  running  average  of  wind. 


31 


RELATIVE  POWER  (DB)  RELATIVE  POWER  < DB )  RELATIVE  POWER  ( DB ) 


1 

I 

k 

| 

! 

i  ! 

'V 

i 

|  i 

^VA..yV^x 

i 

i 

_ _ 

| 

(c)  RANGE  CELL  4 

- 1 1 

2  4 

-  5 

62 

0 

0 

5 

62 

I  l  24 

FREQUENCY  (HZ) 
SPECTRAL  ESTIMATION 


FREQUENCY  (HZ) 
SFECTRAL  ESTIMATION 


Fig.  IJ-20.  Clutter  from  area  2.  OFT  length  =  256,  resolution  =  0.165  Hz,  sampling 
rate  -  28.12  Hz.  (a)  DFTs  averaged  =  265;  (b)  through  (f)  DFTs  averaged  =  53. 
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Fig.  11-21.  Clutter  from  area  2.  DFT  length  -  256,  resolution  0. 165  ITz,  sampling 
rate  -  28.12  Hz.  (a)  DFTs  averaged  =  265;  (b)  through  (f)  DFTs  averaged  -  5  3. 
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Fig.  11-22.  Clutter  from  area  2;  wind  speed  at  antenna  tower,  (a)  Running  standard 
deviation  of  wind;  (b)  running  average  of  wind. 
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Fig,  11-23.  Clutter  from  area  3.  OFT  length  256,  resolution  0. 0412  Hz,  sampling 
rate  =  7.03  Hz.  (a)  DFTs  averaged  375;  (b)  through  (f)  DFTs  averaged  =  75. 
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Fig.  11-24.  Clutter  from  area  3.  DFT  length  256,  resolution  =  0.  0412  Hz,  sampling 
rate  =  7.03  Hz.  (a)  DFTs  averaged  =  375;  (b)  through  (f)  DFTs  averaged  =  75. 
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Fig.  II-25(a-e).  Clutter  from  area  3;  correlation  between  range  cells. 
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Fig.  11-26.  Clutter  from  area  3;  wind  speed  at  antenna  tower,  (a)  Running  standard 
deviation  of  wind;  (b)  running  average  of  wind. 
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Fig.  11-27.  Clutter  from  area  3.  DFT  length  256,  resolution  -  0.0659  Hz,  sampling 
rate  =  11.25  Hz.  (a)  DFTs  averaged  =  540;  (b)  through  (f)  DFTs  averaged  =  108. 
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Fig.  11-28.  Clutter  from  area  3.  DFT  length  =  256,  resolution  =  0.0659  Hz,  sampling 
rate  =  11.25  Hz.  (a)  DFTs  averaged  =  540;  (b)  through  (f)  DFTs  averaged  =  108. 
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Fig.  II-29(a-e).  Clutter  from  area  3;  correlation  between  range  cells. 
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Ill.  TARGET  SPECTRUM  FROM  FOLIATED  AREAS 


R.  D.  Yates 


A.  Theory 

In  a  forest,  the  incident  electromagnetic  energy  is  scattered  from  many  points  which  are 
moving  in  a  bound  manner  due  to  the  wind  turbulence.  This  produces  an  electromagnetic  field 
across  a  range-azimuth  cell  which  is  both  spatially  and  time  varying.  As  a  target  travels 
through  the  range  cell,  this  varying  electromagnetic  field  modulates  its  return  and  produces  an 
observed  signal  with  a  bandwidth  greater  than  the  bandwidth  of  the  same  target  return  in  a  clear 
region. 

A  detailed  study  of  this  problem  is  presented  in  Ref.  6.  The  purpose  of  this  section  is  to  re¬ 
cast  the  results  in  radar  terminology  and  to  summarize  them.  The  forest  model  used  is  identical 
to  that  used  to  model  the  clutter  returns  as  described  in  Sec.  II- A. 

A  point  target  with  a  polarizability  factor  J’  is  assumed  to  be  traveling  in  a  radial  direction 
toward  the  radar  with  a  constant  velocity  v.  Although  a  general  expression  for  the  target  spec¬ 
trum  is  derived,  it  will  be  assumed  that  the  random  permittivity  £(r,  t)  is  statistically  homo¬ 
geneous  and  isotropic  in  all  directions,  and  statistically  stationary  in  time. 

Let  X(t)  =  I(t)  +  iQ( t )  represent  the  complex  envelope  of  the  range-gated  video  signal  from  a 
pulse-doppler  radar  such  as  the  Long  Range  Demonstration  Radar  (LRDR).  Assuming  that  the 
compressed  pulse  duration  is  greater  than  the  time  spreading  due  to  multipath,  the  return  from 
the  point  target  may  be  expressed  as 

X(t)=Gr1/2-  P^/Z  ■  Pf^(t-t0»-  a(t)  (III-l) 


in  which 

=  receiver  gain 
=  target  received  peak  power 

P(t)  compressed  pulse  which  defines  the  range  cell  with 
max  |  P(t)  |  =  1 

a(t)  unit  power  modulation  process  due  to  target  motion  and 
varying  field  across  range  cell. 

In  general,  the  random  process  a(t)  will  not  be  statistically  stationary.  However,  if  the 
random  permittivity  c(r,  t)  is  statistically  homogeneous  and  isotropic  in  r_  and  statistically  sta¬ 
tionary  in  t,  then  as  long  as  the  target's  velocity  is  constant,  a(t)  will  be  a  stationary  random 
process.  For  a  target  traveling  at  a  constant  velocity  v  in  a  radial  toward  the  radar,  in  a  clear 
area,  a(t)  is  then  equal  to  exp  {j  [(4tt/A)  vt  +  0]},  where  0  is  a  random  variable  uniformly  dis¬ 
tributed  between— 7T  and  tt  radians. 

The  function  P  [(2v/c)  t]  represents  a  time  gating  pulse,  and  corresponds  to  the  complex 
envelope  of  the  compressed  transmitted  signal.  In  the  LRDR,  since  the  time  sidelobes  of  the 
compressed  pulse  are  at  least  30  dB  below  the  main  peak,  P(t)  can  be  taken  as  the  complex  enve¬ 
lope  of  the  signal  that  is  fed  to  the  mixers  when  the  input  signal  is  the  basic  pulse  that  defines  a 
bit  in  the  code. 

It  has  been  assumed  that  the  time  multipath  spreading  of  the  basic  bit  pulse  is  small.  If 
this  is  not  the  case,  then  the  duration  of  X(t)  will  be  greater  than  P  [(2v/c)  t].  However,  for  a 
homogeneous,  isotropic  forest,  the  energy  spectral  density  of  X(t)  will  not  be  changed. 
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( III—  2 ) 


The  targets'  received  peak  power  Prp  can  be  expressed  as 
Prp  =  *  G  •  A  •  (Trp  F,p 


in  which 


P^  =  transmitted  peak  power 


G  =  one-way  antenna  gain 

2 

A  =  antenna  aperture  -  GX  /4n 


o’ rp  -  target  cross  section  in  square  meters 


Frp  =  two-way  target  propagation  factor. 


The  antenna  gain  G  has  been  defined  in  Sec.  II- A,  Eq.  (II-5).  As  for  the  clutter  case,  aTk  T  can¬ 
not  be  factored;  their  product  is  given  as 


(HI-  3) 


in  which 


h^  =  antenna  height 
hj.  -  forest  height 

A  =  free- space  wavelength 
R  -  range  to  cell 

=  free- space  permittivity 

■q  =  slab  refractive  index 
g(h)  =  ground  lobing  and  attenuation  factor 
hT  =  target  height 

Pag  =  power  in  the  random  component  of  the  target  return  due 
to  multipath 

f  -  targets'  polarizability. 


The  function  g(h)  represents  the  effects  of  attenuation  and  ground  lobing  within  the  dielectric 
slab  and  is  given  by  Eq.  (II-9).  Although  it  is  not  possible  to  factor  the  product  F  Eq.  (Ill- 3) 
contains  three  propagation  terms  as  given  by  Eq.  (11-17). 

The  function  a(t)  in  Eq.  ( III- 1 )  represents  the  modulation  due  to  the  target's  motion  and  vary¬ 
ing  field  across  the  range  cell,  and  has  a  power  spectral  density  8^(0  as  follows: 


6(f-V+SAC<f> 


(III- 4 ) 
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in  which 

=  dopplcr  frequency  =  2v/A 
6(f)  =  delta  function 

SAc<f)  =  multipath  contribution  spectrum 

pAc=r  sAC<f)rff. 

^  -  oo 


In  Eq.  ( III— 4 ),  the  spectrum  of  the  direct  path  to  the  target  is  the  impulse  at  the  doppler  fre¬ 
quency  Ej,  while  the  contribution  due  to  the  multipath  is  given  by  S^^(f).  The  spectrum  S^(f)  of 
the  random  process  a(t)  has  been  normalized  to  unit  power. 

The  multipath  component  spectrum  can  be  written  as 


sAC<f) 


■£ 


SQ(u)  Sw(  ^  r?  n/1  -  u,  f  -  rjfdu)  du 


( III— 5 ) 


in  which 

s  (u)  =  -  •  PN i  -  u)  •  <p(  v  'ft  -  u ) 

X  €  a 
o 

w’  here 

p{x)  ~  ^  (3  +  2x^  +  3x^) 

and  <p(k)  is  the  spatial  spectrum  of  the  random  permittivity  and  is  defined  as 


(III-6) 


( HI-  7 ) 


C(r)  e — d 


3 

r 


=  (p(  |k|)  due  to  statistical  homogeneity  and  isotropy  (III- 8 ) 

where  C(r)  is  the  correlation  function  of  the  random  permittivity  and  is  defined  as 

C(r)  -  E{*(R,t)  e*(R-  r,  t)} 

=  C(  Ir  |)  due  to  statistical  homogeneity  and  isotropy  and 

temporal  stationarity  .  (III-9) 

The  function  ( k,  f)  in  Eq.  (Ill- 5)  represents  the  effects  of  the  motion  of  the  scatterers  due 
to  the  velocity  field  V(r,  t),  and  can  be  expressed  as 

Sw(k,  f)  =  f  q(k,  t)  e~iZ*tT  dT  .  (III-10) 

in  which 

q(k,  r)  =  E(exp{jk  [D(t)-  D(t  -  r)]})  .  (III-ll) 

In  Eq.  (Ill— 1 1 ),  D(t)  is  a  displacement  random  process  and  is  not  a  function  of  position,  since  the 
forest  has  been  taken  to  be  statistically  homogeneous  and  isotropic  in  all  directions.  The  cor¬ 
relation  function  m(r)  defined  by  Eq.  (11-28)  is  seen  to  be  equal  to  q  [(4tt /\),t].  Hence,  the  unit 
power  clutter  spectrum  is  then  equal  to  [(47r/A),  f]. 
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The  first  simplification  of  Eq.  (Ill- 5)  is  when  the  scatterers  are  not  in  motion,  so  that 
Sw(k,  f)  -  6(f).  Hence,  the  power  spectral  density  SA(,(f)  given  by  Eq.  (ill- 5)  reduces  to 


SAC(f) 


for  0  <  f  ^  r7  f ^ 


=  0  elsewhere 

Thus,  the  power  spectral  density  of  the  random  process  a(t)  can  he  written  as 


(III-  1  2 ) 


sa(f)  = 


1  +  p 


AC 


6(f  -  fd)  + 


*  ft 


for  0  <  f  <  ??fd 


(III-  1  3) 


in  which  |3(x)  and  c?(k)  are  defined  by  Eqs.  ( III  -  7 )  and  ( III-  8 ),  respectively. 

It  is  seen  that  the  target  spectrum  contains  no  frequencies  above  the  doppler  frequency 
because  the  target  is  traveling  in  a  radial  direction  toward  the  radar.  The  power  due  to 

the  multipath  is  inversely  proportional  to  a,  the  attenuation  factor  within  the  forest.  Hence,  for 
a  forest  with  thick  underbrush,  one  would  expect  the  impulse  component  in  Eq.  ( III- 1  3 )  to  dominate. 
For  a  telephone-type  forest  with  very  little  underbrush,  the  random  component  S^p(f)  would  be 
the  dominant  term. 

Since  the  observed  signal  X(t)  given  by  Eq.  ( I II - 1 )  is  a  time-gated  version  of  the  process  a(t), 
we  must  consider  the  energy  spectral  density  of  X(t)  defined  as  follows: 


Ex(f)  =  E  { I  \  x<t)  e 


r 


jzxft  dt 

-  b  •  \  SaO>)|F(f-  v)\2 


( III-  1 4 ) 

(III-  1  5 ) 


where  b  is  a  constant,  and  F(f)  is  the  Fourier  transform  of  the  gating  pulse  used.  Hence,  the 
resolution  of  the  random  component  S^^(f)  of  a(t)  will  depend  on  the  ratio  of  the  bandwidth  of 
SAr(f)  to  the  bandwidth  of  |F(f)|^.  If  this  ratio  is  large,  then  Ex(f)  can  be  approximated  as 
follows: 


Ex(f)  -b'  [W(f-fd>  +  SAC(f)] 

in  which  b'  is  a  constant,  and  W(f)  is  a  unit  area  window  function  defined  as 

w(f)  =  — [nnl2— 


(III- 16 ) 


(III-  17 ) 


/  iFlf )  |  df 

The  other  simplification  of  Eq.  (Ill- 5)  results  when  the  target  velocity  is  assumed  to  be  zero 
and  wind  turbulence  is  present.  Hence, 


SAC(f)  =  )  So(u)  SW(  T  V  ^  -  u>  f)  du 


(111-18) 


The  normalized  clutter  power  spectral  density  is  given  as  Sw  [ ( 4 7r/A ),  f ] .  Assume  that  the 
phase  modulation  index  \±  ~  Ana/X,  given  by  Eq.  (11-20),  is  less  than  1.  Then,  Eq.  (111-4)  can  be 
expressed  as 
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( III- 1 9  ) 


Sa(f>  (1  +  PACrl  '  II1  +  PAc(i  "  ^)|  6(f)  +  PAC  '  ^Sd(f)| 

in  which  6(f)  is  the  delta  function,  S , (f )  is  the  normalized  displacement  process  defined  by 
“2  d 

Eq.  (11-40),  and  p  is  a  mean  square  phase  modulation  index  defined  as  follows: 


2 

P 


47T(7  n  ^2 

X  } 


/q1  uSo(1  ~  U)  du 

fo  So(u)du 


( III—  2  0) 


The  function  S^(u)  is  given  by  Eq.  (III-6),  and  [the  power  in  the  AC  component  of  a ( t ) )  is 

given  as 


PAC-I  So(u)du 


( II I- 2  1 ) 


2  2 

It  follows  from  Eq.  (III-20)  that  p  <p  ,  where  p  is  the  phase  modulation  associated  with  the 
clutter  spectrum.  Since  the  non-impulse  part  of  Eq.  (Ill— 1 9 )  is  due  mainly  to  forward  scattering, 
one  would  expect  that  the  DC-to-AC  power  ratio  in  S^(f)  would  be  greater  than  that  of  the  clutter 
spectrum. 


B.  Experiments 

In  this  section,  several  experiments  are  described  which  were  performed  to  verify  the 
theoretical  results  given  in  Sec.  A  above. 

1.  Transponder  Returns 

A  coherent  transponder  described  in  Sec.  Ill- F- 1  in  Vol.  1  was  used  as  a  point  target.  Its 
antenna  is  horizontally  polarized  and  essentially  omnidirectional.  The  received  signal  is  ampli¬ 
fied,  delayed  by  20psec,  and  frequency  shifted  by  (256  x  555.6  x  10  *  7.03  Hz.  The  reason 

for  the  frequency  shift  is  to  ensure  that  the  target  return  does  not  overlap  the  clutter  return  in 
frequency.  Amplification  of  the  signal  was  necessary  in  order  to  ensure  that  the  target  spectrum 
be  well  above  the  receiver  noise.  The  transponder  was  pulled  at  a  constant  velocity  at  a  fixed 
height  above  the  nominally  flat  ground  within  foliated  areas. 

a.  Transponder  Closed- Loop  Stability 

In  order  to  determine  the  spectral  purity  of  the  transponder,  a  closed-loop  stability  test  was 
performed.  The  low  level  coded  gated  RF  from  the  LRDR  was  fed  directly  to  the  transponder 
via  a  circulator.  The  amplified,  time,  and  frequency- shifted  signals  were  then  fed  to  the  LRDR 
receiver  input.  Adequate  attenuation  was  used  to  ensure  that  the  transponder  and  LRDR  receiver 
were  operating  linearly.  Using  the  measurement  mode  of  the  LRDR,  the  compressed,  range¬ 
gated,  low-pass-filtered,  quadrature  video  components  were  recorded  on  digital  magnetic  tape 
using  the  TAPWRITE  system  described  in  Sec.  III-F-3  in  Vol.  1. 

An  estimate  of  the  power  spectral  density  of  the  transponder  signal  is  given  in  Fig.  Ill- 1 . 

The  main  line  is  at  +7.03  Hz  while  its  image  at  —7.03  Hz  is  28  dB  below  it.  This  image  is  due  to 
phase  and  amplitude  errors  in  the  quadrature  channels.  The  other  lines,  which  are  at  least 
38  dB  below  the  main  line,  appear  to  be  due  to  the  transponder’s  nonideal  digital  frequency  shifter. 
A  closed-loop- stability  test  of  the  entire  radar,  as  reported  in  Sec.  II- D  in  Vol.  1,  shows  that 
there  are  no  spurious  signals  present.  These  extra  lines  presented  no  major  problems  in  the 
experiments  to  be  discussed. 
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b.  Spectrum  Broadening  Due  to  Multipath 

In  order  to  verify  the  nature  of  the  power  spectral  density  [as  given  by  Eq.  { III- 1  3 ) ]  for  the 
return  from  a  point  target  in  a  foliated  area,  the  transponder  was  used  as  the  target. 

The  area  chosen  was  several  kilometers  from  the  radar  and  consisted  of  trees  with  an  aver¬ 
age  height  of  40  feet  with  light  foliage.  The  transponder  traveled  a  distance  of  9  0  feet  in  which 
a  50-foot  range  cell  was  centered.  Eleven  round  trips  were  made  so  that  some  averaging  of  the 
data  was  possible. 

The  time  data  per  pass  were  examined  to  isolate  the  section  which  contained  the  major  re¬ 
turn.  A  time-weighted  FFT  was  performed  on  the  data  and  successive  FFTs  were  added  inco¬ 
herently  to  provide  smoothing.  If  the  number  of  points  per  FFT  is  denoted  as  N  -  2n,  then  the 
number  of  FFTs  per  pass  was  1024/N. 

The  estimated  spectrum  for  an  inbound  pass  is  shown  in  Figs.  Ill- 2,  -3,  and  -4  for  FFT 
lengths  of  64,  256,  and  1024,  respectively.  The  corresponding  spectrum  for  an  outbound  pass 
is  shown  in  Figs.  Ill- 5,  -6,  and  -7.  It  appears  that  Figs.  Ill- 3  and  -6  show  the  best  tradeoff  be¬ 
tween  resolution  and  smoothing  and,  on  examining  them,  it  is  seen  that  the  inbound  and  outbound 
spectra  are  mirror  images  of  each  other  about  the  offset  frequency  of  7.03IIz.  The  predicted 
asymmetry  about  the  offset  doppler  frequency  of  8.21  IIz  for  the  inbound  pass  and  5.85  Hz  for  the 
outbound  pass  is  also  present.  Since  the  target  doppler  is  1.18  Hz  and  the  3-dB  bandwidth  about 
1.4  Hz,  the  Q  of  the  spectrum  is  then  0.85.  Although  this  appears  to  be  smaller  than  expected, 
it  does  correlate  with  the  fact  that  the  detection  of  a  single  man  in  the  same  range  cell  was  very 
difficult.  It  appears  that  the  multipath  was  severe  enough  to  cause  considerable  broadening  of 
the  target  spectrum. 

In  all  these  Figures,  the  component  about  zero  frequency  is  the  spectrum  of  the  clutter  from 
the  range  cell  20psec  further  out.  The  component  about  —7  Hz  is  due  to  the  sideband  unbalance 
of  the  system.  The  wind  in  the  area  of  the  range  cell  was  minimal  so  that  spectral  broadening 
due  to  the  motion  of  the  scatterers  can  be  neglected. 

c.  Spectrum  Broadening  Due  to  Wind 

In  order  to  determine  the  effect  of  the  motion  of  the  scattering  centers  on  the  spectra  pre¬ 
sented  in  Sec.  b  above,  the  transponder  was  placed  in  the  center  of  the  range  cell  and  the  result¬ 
ing  radar  return  was  recorded.  The  theoretically  predicted  spectrum  is  given  by  Eq.  (Ill- 19). 

The  estimated  spectrum  is  given  in  Fig.  Ill- 8,  while  the  running  mean  and  standard  deviation 
of  the  wind  velocity  during  the  data  recording  time  are  given  in  Fig.  III-9(a-b).  Wind  velocity 
was  measured  by  an  anemometer  at  the  top  of  the  antenna  tower,  and  gives  a  relative  measure 
of  the  wind  conditions  at  the  range  cell  since  it  was  several  kilometers  away. 

It  is  seen  that  the  spectral  broadening  about  the  transponder  offset  frequency  of  7.03  IIz  occurs 
about  25  dB  below  the  main  line.  Except  for  the  clutter  spectrum  about  zero  frequency,  the  other 
lines  present  are  due  to  the  transponders  digital  frequency  shifter  as  discussed  in  Sec.  a  above. 

Based  on  the  above  results,  it  is  seen  that  the  effect  of  the  wind  on  the  spectra  presented  in 
Sec.b  can  be  ignored.  However,  for  high  wind  velocities,  this  may  not  be  the  case  since  the 
DC-to-AC  power  ratio  of  the  spectrum  given  in  Fig.  Ill—  8  would  be  greatly  reduced. 

2.  Spectrum  of  Man  in  Woods 

A  man  is  an  extended  target  vertically,  so  that  some  form  of  averaging  of  the  field  will  occur. 
If  the  field  statistics  do  not  change  significantly  over  the  target's  height,  the  spectrum  of  a  man 
should  be  similar  to  that  of  a  point  target  such  as  a  transponder. 
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One  man  walked  at  a  constant  velocity  in  a  radial  direction  with  respect  to  the  radar.  The 
area  consisted  of  trees  with  light  underbrush.  Several  passes  were  made  through  the  range  cell 
to  ensure  smoothing  of  the  spectral  estimates.  The  range  cell  was  several  kilometers  from  the 
radar  and  in  a  different  location  from  the  transponder  test  area. 

The  spectral  estimates  with  a  frequency  resolution  of  0.659,  0.11,  and  0.0549  Hz  are  given 
in  Figs.  111-10,  -11,  and  -12  for  an  inbound  pass,  and  in  Figs.  HI-13,  -14,  and  -15  for  an  outbound 
pass.  It  is  clear  that  a  frequency  resolution  of  0.659  Hz  is  not  adequate  in  order  to  estimate  the 
target  spectrum.  However,  this  case  is  presented  here  for  comparison  since  the  bandwidth  of 
the  coherent  filters  in  the  LRDR  signal-processing  scheme  is  0.659  Hz.  It  appears  that  the  best 
estimates  are  given  by  Figs.  111-12  and  -15. 

On  examining  Fig.  Ill- 12  for  the  inbound  pass,  the  definite  asymmetry  about  the  target  doppler 
of  4.  3  Hz  is  present.  The  3-,  6-,  and  1  0-dB  bandwidth  measures  are  0.42,  0.93,  and  1.46  Hz, 
respectively.  Hence,  the  corresponding  Q s,  defined  as  the  ratio  of  the  doppler  to  bandwidth, 
are  10,  4.6,  and  3,  respectively.  Since  the  bandwidth  is  proportional  to  the  target’s  velocity, 
the  Q  is  independent  of  target  velocity. 

The  outbound  spectrum  is  not  a  mirror  image  of  the  inbound  spectrum  as  predicted  by  the 
theory  and  the  transponder  experiments.  Due  to  the  limited  time  allowed  for  science  experiments, 
this  problem  was  not  investigated  further.  It  could  be  attributed  to  poor  velocity  control  of  the 
target,  and  an  incorrect  centering  of  the  range  cell  and  data  window. 

3.  Recommendations  for  Further  Work 

An  extensive  study  of  the  spectrum  of  a  target  in  a  clear  region  was  presented  in  Ref.  5.  In 

addition,  the  spectrum  of  a  man  in  foliage  from  one  range  cell  was  presented  and  it  showed  the 

-  8 

asymmetry  about  the  doppler  frequency.  The  R  dependence  of  the  received  target  power  on 
range  from  foliated  areas,  as  given  by  Eq.  (II- 3),  was  also  verified.  The  refractive  index  r?  and 
attenuation  cv  for  several  types  of  foliage  were  also  estimated. 

In  this  report,  the  transponder  returns  appear  to  have  a  spectrum  with  the  correct  properties 
as  predicted  by  Eq.  (11-13). 

An  extensive  experimental  program  should  be  carried  out  to  check  the  validity  of  the  expres¬ 
sions  for  the  target  spectrum  as  presented  in  Sec.  III-A.  Different  types  of  foliage  should  be 
studied  to  determine  how  they  affect  the  target  spectrum.  Assuming  that,  over  the  observation 
time,  the  target's  velocity  is  constant,  it  should  be  possible  to  characterize  the  foliage  by  a  Q 
which  can  then  be  used  in  the  design  of  the  target  filter  bank. 

Since  the  target  return  consists  of  a  specular  plus  a  random  component,  experiments  should 
be  designed  to  measure  the  ratio  of  the  specular  to  random  power.  If  the  specular  power  is 
dominant,  then  the  receiver  should  coherently  integrate  over  the  time-on-target.  However,  if 
the  random  component  is  dominant,  then  the  receiver  should  coherently  integrate  over  its  coher¬ 
ence  time  and  incoherently  integrate  over  the  time-on-target.  This  assumes  that  the  target’s 
velocity  is  constant  over  the  time-on-target.  Variations  of  target  velocity  will  spread  the  ap¬ 
propriate  spectrum  further,  and  require  wider  bandwidth  filters. 

In  order  to  predict  receiver  performance,  the  statistics  of  the  target  return  must  be  known. 

If  there  is  a  great  deal  of  multipath  present,  one  might  speculate  that  the  random  component  of 
the  target  return  is  a  Gaussian  process.  Experiments  should  be  designed  to  estimate  the  statis¬ 
tics  of  the  target  return  for  various  types  of  foliage. 

Finally,  more  experiments  are  needed  to  determine  the  behavior  of  the  received  target 
power  as  a  function  of  the  type  of  foliage. 
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Fig.  Ill- 1 .  Transponder's  closed- loop- stability  spectrum. 
DFT  length  -  256,  sampling  rate  28.12  Hz,  resolution  - 
0.165  Hz,  DFTs  averaged  =  54. 


Fig.  Ill- 2 .  Spectrum  of  transponder  in  woods  inbound.  DFT  length  -  64, 
sampling  rate  =  28.12  Hz,  resolution  =  0.659  Hz,  DFTs  averaged  -  176, 
target  passes  =11,  target  velocity  =  1.35  ft/sec. 
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Fig.  Ill- 3.  Spectrum  of  transponder  in  woods  inbound.  DFT  length  =  256, 
sampling  rate  =28. 12  Hz,  resolution  =  0.165  Hz,  DFTs  averaged  =  44, 
target  passes  :  11,  target  velocity  =  1.35  ft/sec. 


FREQUENCY  (Hz) 


Fig.  Ill—  4.  Spectrum  of  transponder  in  woods  inbound.  DFT  length  =  1024, 
sampling  rate  =  28.12  Hz,  resolution  =  0.0412  Hz,  DFTs  averaged  =  11, 
target  passes  =  11,  target  velocity  =  1.35  ft/sec. 
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Fig.  Ill- 5.  Spectrum  of  transponder  in  woods  outbound.  DFT  length  64, 
sampling  rate  28.12  Hz,  resolution  0.659  Hz,  DFTs  averaged  176, 
target  passes  11,  target  velocity  -  1.35  ft/s ec. 


Fig.  Ill- 6 .  Spectrum  of  transponder  in  woods  outbound.  DFT  length  256, 
sampling  rate  28.12  Hz,  resolution  =  0.165  Hz,  DFTs  averaged  44, 
target  passes  =11,  target  velocity  =  1.35  ft/s  ec. 
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Fig.  III-7.  Spectrum  of  transponder  in  woods  outbound.  DFT  length  =  1024 
sampling  rate  =  28.12  Hz,  resolution  =  0.0412  Hz,  DFTs  averaged  =  11 
target  passes  =  11,  target  velocity  =  1.35  ft/sec. 


Fig.  Ill- 8.  Spectrum  of  transponder  at  rest  in  woods  with  wind. 
DFT  length  =  256,  sampling  rate  =  28.12Hz,  resolution  =  0.165  Hz, 
DFTs  averaged  =  54. 
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Fig.  Ill— 9 .  Transponder  at  rest  in  woods  with  wind,  (a)  Running  standard 
deviation  of  wind;  (b)  running  average  of  wind. 
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Fig.  III- 10.  Spectrum  of  man  in  woods  inbound.  DFT  length  =  64,  sampling 
rate  -  28.12Hz,  resolution  =  0.659  Hz,  DFTs  averaged  =  132,  target  passes  =  11, 
target  velocity  1.5  m/sec. 


Fig.  IH-11.  Spectrum  of  man  in  woods  inbound.  DFT  length  =  256,  sampling 
rate  18.  75  Hz,  resolution  =0.11  Hz,  DFTs  averaged  -10,  target  passes  =  10, 
target  velocity  =  1.5  m/sec. 
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Fig.  II J- 1 2.  Spectrum  of  man  in  woods  inbound.  OFT  length  512.  sampling 
rate  18.75  Hz,  resolution  0.0549  Hz,  DFTs  averaged  10,  target  passes  10, 
target  velocity  1.5  m/sec. 


Fig.  Ill- 13.  Spectrum  of  man  in  woods  outbound.  DFT  length  64.  sampling 
rate  28.12Hz,  resolution  0.659  Hz,  DFTs  averaged  132,  target  passes  11, 
target  velocity  =  1.5  m/sec. 
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Fig.  Ill- 14.  Spectrum  of  man  in  woods  outbound.  DFT  length  =  256,  sampling 
rate  =  18.75  Hz,  resolution  =  0.11  Hz,  DFTs  averaged  =  10,  target  passes  =  10, 
target  velocity  =  1.5  m/sec. 


Fig.  Ill- 15.  Spectrum  of  man  in  woods  outbound.  DFT  length  -  512,  sampling 
rate  =  18.75  Hz,  resolution  =  0.0549  Hz,  DFTs  averaged  -  10,  target  passes  =  10, 
target  velocity  1.5  m/sec. 
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IV.  CLUTTER  STATISTICS 


J  .11.  Teele 
R.  D.  Vates 


A.  Theory 

In  order  to  design  optimum  receivers  and  analyze  their  performance,  the  clutter  statistics 
must  be  known.  In  Sec.  II-A,  a  clutter  model  was  presented  which  could  predict  the  clutter 
spectrum.  However,  nothing  was  said  about  the  clutter  statistics.  If  the  correlation  distances 
of  the  random  permittivity  £(r,t)  in  the  transverse  directions  are  denoted  as  r  and  r  ,  then  as 

x  y 

long  as  r  •  r  «  R  •  AR  •  AO,  there  w'ill  be  many  independent  scattering  centers  in  the  range - 

x  y 

azimuth  cell  of  area  R  •  AR  .  AO.  Assuming  that  c(r,  t)  is  homogeneous  and  isotropic  in  the 
transverse  directions  so  that  there  are  no  dominant  scatterers,  one  would  then  postulate  that 
the  clutter  return  is  a  Gaussian  random  process. 

Consider  the  clutter  return  from  a  single  range-azimuth  cell  as  given  by  Eq.  (II-l).  It  has 
a  power  spectral  density  with  an  impulse  at  zero  frequency.  Hence,  the  quadrature  components 
contain  a  DC  or  specular  component  and  may  be  expressed  as 

I(t)  =  G*'2  [IQ  +  Ij(t)l  ( IV  - 1 ) 

1/2 

Q(t)  =  Gr  (Q0+Q1(t)l  .  (IV -2) 

The  random  variables  I  and  represent  the  DC  clutter  component  and  are  zero  mean  and 
uncorrelated  with 

I  -  A  cose  (IV-  3) 

Q  -  A  sin  0  ( IV-4) 

o 

where  0  is  a  random  variable,  independent  of  A,  and  uniformly  distributed  between  0  and  2  jr. 

The  mean  square  value  of  the  random  variable  A  is  the  DC  clutter  power  and  is  given  as 

PDC  =  E<A2}  =  E{lo2  +  Q2} 

=  pc‘Rc(00)  •  (IV-5) 

The  random  processes  I ^ (t)  and  Q ^ (t )  arc  zero  mean,  uncorrelated,  wide-sensc  stationary 
with  correlation  functions 

Rj  (t)=Rq  (r)a|-  Pc-  [Rc(t)-Rc(«)J  .  (IV -6) 

1  *  1 

The  AC  clutter  power  is  then  given  as 

pAC  =E{l2(t)  +  Q2(t)}  =  Pc-PDC  •  (IV -7) 

The  total  clutter  power  is  given  by  Eq.  (11-2),  while  the  normalized  clutter  correlation  func¬ 
tion  Rc<r)  is  given  by  Eq.  (11-15).  In  addition,  the  random  variables  I  and  are  uncorrelated 
with  the  random  processes  1^  (t )  and  Q^(t). 

If  the  clutter  can  be  modeled  as  a  Gaussian  random  process,  then  I  and  Q  are  independent 
Gaussian  random  variables,  A  is  a  Rayleigh  random  variable,  and  I^(t)  and  Q^(t)  are  independent 
Gaussian  random  processes. 
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It  should  be  noted  that  I  and  Q ^  are  spatial  random  variables  so  that  the  expectation  in 
Eq.  ( I V - 5 )  is  with  respect  to  the  ensemble  of  possible  range  cells.  The  random  processes  I^(t) 
and  Q^(t)  evolve  in  time  so  that  the  expectation  in  Eq.  (IV -6)  is  a  temporal  one. 

Now  there  are  two  extensions  of  the  above  model  which  occur  in  practice.  It  was  assumed 
that  the  velocity  field  was  a  stationary  process  in  time.  A  modulation  index  was  defined  in 
Eq.  (11-44)  as 

n(z)  =  ^  •  /S(vo>  •  Z2  (IV-8) 


where  0(v^)  is  a  monotonically  increasing  function  of  the  average  wind  velocity  vq  above  the 
trees.  Let  v  be  slowly  varying  compared  with  the  AC  clutter  signals  and  let  it  have  a  correla¬ 
tion  time  Tq.  Then,  over  a  time  span  T  ,  the  clutter  processes  I(t)  and  Q(t)  can  be  modeled 
as  stationary  and  Gaussian,  conditioned  on  vq  being  known.  Over  periods  much  longer  than  T^, 
the  clutter  probably  is  not  Gaussian  and  certainly  is  not  statistically  stationary. 

The  second  extension  concerns  the  spatial  statistics  of  the  clutter.  The  forest  was  assumed 
to  be  homogeneous  and  isotropic  in  the  transverse  directions.  This  implies  that  the  clutter  area 
cross  section  is  constant.  Suppose  that  cr  is  a  random  process  with  respect  to  the  trans¬ 
verse  directions  x,  y,  and  that  over  the  range  cell  area  R  *  AR  ♦  A0  it  can  be  considered  a 
constant.  Then,  conditioned  on  (tq  being  known,  the  clutter  is  a  Gaussian  process.  However, 
over  the  ensemble  of  possible  range  cells,  the  clutter  is  not  a  Gaussian  process  and  its  statistics 
will  depend  on  the  statistics  of  ctq. 

An  additional  complication  is  that  the  two-way  clutter  propagation  factor  Fc  may  also  he  a 
function  of  space  due  to  different  frce-space  forest  interface  losses,  hills,  etc. 

The  spatial  statistics  of  the  clutter  are  important  because  they  indicate  the  dynamic  range 
needed  in  the  receiver  front  end. 

The  type  of  processing  of  the  time  data  in  a  given  range-azimuth  cell  will  depend  on  the 
temporal  clutter  statistics.  Consider  the  case  of  a  small  phase  modulation  index.  Then,  from 
Eq.  (11-22)  the  clutter  power  and  clutter  correlation  functions  may  be  written  as 


R.  (T)  =  R  (T)  -  4 

h  2 


AC 


pd{t) 


(IV-10) 


AC  " 


where 


2 


is  a  mean  square  modulation  index  given  by  the  expression 


—  ,hf 

n2  (4^)2  1  W(z)  Z4  dz  •  02(vf)) 


(IV-11) 


(IV-12) 


and  Pj)(T)  is  the  normalized  displacement  correlation  function  given  by  Eq.  (11-39). 

It  is  seen  that  the  total  clutter  power  does  not  depend  on  the  average  wind  velocity  v  .  How- 

2  3  ° 

ever,  in  Sec.  II-B  it  was  postulated  that  p  (vq)  was  proportional  to  v  .  Hence,  the  AC  clutter 

power  is  very  sensitive  to  the  wind  velocity.  If  the  Gaussian  nature  of  the  clutter  is  to  be  deduced 

from  studying  the  time  signal  from  one  range  cell,  it  is  important  to  properly  normalize  the 

clutter  so  that  it  has  a  constant  mean  and  variance. 
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Note  that  Eqs.  (IV-9),  (IV-10),  and  (IV-11)  imply  that  the  clutter  spectrum  shape  does  not 
depend  on  the  average  wind  speed  v  ,  only  the  powers  do. 

13.  Experiments 

The  experiments  performed  to  verify  the  clutter  statistics  were  limited  due  to  time.  No 
spatial  statistics  were  studied.  However,  in  Ref.  5  it  was  shown  that,  for  flat,  uniform-type 

-7 

foliage,  the  R  dependence  of  the  clutter  power  is  correct.  The  clutter  return  as  a  function 

of  range,  under  the  absence  of  wind,  was  recorded.  A  least-squares  fit  of  an  R“n  dependence 
Z  Z 

of  I  (R)  +  Q  (II)  was  found.  Due  to  recording  limitations,  the  clutter  statistics  were  not  studied. 

Using  the  measurement  mode  of  the  ERDR,  the  quadrature  components  of  the  clutter  from 
one  range  cell  were  recorded  for  an  hour.  One  quadrature  channel  was  studied  and  the  data 
reduced  to  several  segments  of  8192  points,  in  each  segment,  the  sampling  rate  was  11.25  Hz 
corresponding  to  a  segment  duration  of  723  sec.  A  histogram  of  one  of  these  segments  of  8192 
points  is  shown  in  Fig.  IV- 1,  where  it  is  seen  that  the  probability  density  is  not  Gaussian.  Similar 
results  were  found  for  other  segments. 

The  data  were  now  filtered  to  remove  the  high-energy  clutter  about  zero  frequency  and  the 

receiver  noise.  A  histogram  of  the  filtered  data  showed  no  improvement.  For  each  segment, 

the  8192  points  were  divided  into  B  consecutive  blocks,  with  8192/B  points  per  block.  Given 

2 

a  block,  the  sample  mean  M  and  sample  variance  s  were  computed  and  used  to  normalize  the 
data  in  the  block  as  follows: 

x.  -  M 

x. (IV  - 13) 

l  s 

where  are  the  original  data  points,  and  x^  are  the  normalized  data.  Histograms  were  now 
made  for  each  segment.  The  case  for  B  32  is  shown  in  big.  IV-2  for  the  same  data  used  in 
Fig.  IV-1.  Each  block  consisted  of  256  points  corresponding  to  a  time  span  of  23  sec.  By  using 

7 

the  Kalmogorov-Smirnov  nonparametric  tests,  the  probability  that  the  8192  points  were  independ- 

2 

ent  Gaussian  with  variance  o  and  mean  m  was  50  oereent. 

Other  block  lengths  were  tried  since  there  is  a  tradeoff  between  the  number  of  samples 
needed  to  estimate  the  mean  and  variance,  and  the  time  over  which  the  clutter  statistics  are 
stationary.  It  was  found  that  B  32  gave  the  highest  probability  of  being  Gaussian. 

If  the  N  -  8192/B  clutter  samples  in  a  block  were  independent  identically  distributed 
Gaussian  random  variables,  then  xL  as  given  by  Eq.  (IV- 13)  would  actually  be  distributed  as 
\l  1  —  1  /N  •  t^  ^ ,  where  t^  has  a  student-t  distribution  with  N  —  1  degrees  of  freedom.  How¬ 
ever,  for  N  >  40  there  is  negligible  difference  between  the  t-density  and  the  Gaussian  density 
over  the  ±3<r  points. 

The  fact  that  the  histogram  in  Fig.  IV-2  has  a  greater  variance  than  one  obtained  from  a 
set  of  8192  identically  distributed,  independent,  Gaussian  random  variables  can  be  caused  by 
the  following  conditions.  The  first  condition  is  that  the  data  points  are  not  independent,  so  that  the 
equivalent  degrees  of  freedom  are  less  than  8192;  the  second  is  due  to  the  nonstationary  clutter 
statistics  so  that  the  data  are  not  normalized  properly. 

Based  on  the  above  results,  it  appears  that  the  first-order  statistics  of  the  clutter  are 
Gaussian  over  the  correlation  time  of  the  average  wind  velocity  above  the  trees. 
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However,  an  extensive  study  should  be  made  of  the  temporal  statistics  of  the  clutter  from 
several  different  range  cells  under  different  wind  conditions.  The  average  wind  velocity  above 
the  trees  in  the  area  studied  should  be  recorded.  This  will  help  in  correlating  wind  changes 
with  those  in  the  clutter  spectral  shape  and  power,  and  in  properly  normalizing  the  data. 

The  spatial  statistics  of  the  clutter  from  foliated  areas  should  also  be  studied.  Let 
I(R,  0,  t  )  and  Q(R,  0,  t  )  be  the  quadrature  components  of  the  clutter  from  a  range-azimuth  cell 
at  R,  O.  The  quantity 


Z(R,  0) 


Y  jTT  ii2(r,  e,  t)  +  q2(r,  e,  t)J  dt 

P  •  G  •  A  •  R  •  AR  •  AO 

O 


(IV- 1 4) 


could  be  computed  where  T  is  long  enough  to  smooth  out  the  variations  due  to  wind.  The  quan¬ 
tities  in  the  denominator  are  defined  in  Eqs.  (II- 2)  and  ( II - 3) . 

If  the  forest  were  homogeneous  and  isotropic  over  the  range-azimuth  cells  measured,  then 
{D(R^,  0^)}  would  be  exponentially  distributed  random  variables  with  a  mean  value  of  cr^F^,  as 
given  by  Eq.  (II- 7).  Hcncc,  a  histogram  of  {L(R^,  0^)  •  R^}  should  be  made  to  see  if  it  is 
exponential. 


Fig.  IV- i.  Histogram  of  unnormalized  clutter  data. 
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Fig.  IV-2.  Histogram  of  normalized  clutter  data. 
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V.  SUBCLUTTER  VISIBILITY  IN  DIGITAL  MTI  RADARS 


L.  Cartledge 
R.  D.  Yates 


The  purpose  of  this  section  is  to  present  an  expression  for  the  subclutter  visibility  in  Digital 
Moving  Target  Indicating  Radars  such  as  the  LRDR. 

Consider  a  range-gated,  pulse-doppler,  MTI  radar  whose  general  structure  is  described  by 
Figs.  V-t  and  V-2.  The  transmitted  signal  consists  of  an  RF  pulse  which  is  sent  at  a  fixed  prf 
(pulse  repetition  frequency).  Generally,  this  pulse  will  be  modulated  so  that  pulse  compression 
can  be  used  to  increase  the  range  resolution.  For  this  discussion,  it  will  be  assumed  that  the 
RF  pulse  consists  of  N  coherent  basic  subpulses.  A  binary  code  is  modulated  onto  the  RF  pulse 
by  letting  the  basic  subpulses  have  an  amplitude  of  ±A.  In  this  manner,  the  range  cell  is  defined 
by  the  basic  subpulse  provided  that  the  code  has  a  correlation  function  with  low  time  sidelobes. 

The  analog  section  of  the  receiver  is  shown  in  Fig.  V-l.  The  first  bandpass  filter  is  used 
to  limit  the  out-of-band  interference.  It  is  followed  by  a  low-noise,  wideband,  RF  amplifier. 

This  amplifier  then  feeds  a  variable  attenuator,  which  is  followed  by  a  second  wideband  low-noise 
RF  amplifier,  a  second  bandpass  filter,  quadrature  video  detectors,  and  video  amplifiers.  In 
practice,  the  combination  of  the  second  RF  amplifier  and  bandpass  filter  is  implemented  with 
several  cascaded  wideband  RF  amplifiers,  the  bandpass  filter  being  placed  in  the  center  of  the 
chain.  In  a  properly  designed  system,  the  video  noise  level  is  set  by  the  second  RF  amplifier 
for  most  of  the  attenuator's  range. 

The  digital  section  of  the  receiver  is  shown  in  Fig.  V-2.  The  in-phase  and  quadrature  video 
signals  are  sampled  by  two  analog-to-digital  (A/D)  converters,  followed  by  a  digital  transversal 
filter  which  performs  the  pulse  compression.  The  compressed  samples  are  then  fed  to  a  digital 
signal  processor  which  produces  a  detection  statistic  for  each  range-doppler  cell  to  be  examined. 
These  statistics  then  are  sent  to  a  thresholding  device  for  decision-making  and  a  display  unit. 

The  digital  signal  processor  may  be  implemented  in  several  ways.  The  most  flexible  is  to 
use  a  general-purpose,  high-speed  digital  computer  which  is  programmed  to  compute  the  detec¬ 
tion  statistics  given  the  input  data.  Signal  processing  algorithms  then  exist  as  software  packages 
which  can  be  easily  modified.  The  most  inflexible  method  is  to  use  a  hardwired  special-purpose 
computer.  Rased  on  the  experience  gained  with  the  LRDR,  the  first  method  should  always  be 
used  until  all  the  details  and  operating  characteristics  of  the  signal  processor  are  understood. 

The  fixed-point  arithmetic  and  finite  word  length  of  the  computing  machine  used  can  produce 
quantization  and  dynamic  range  problems  which  can  affect  the  performance  of  the  signal  processor 
in  unpredictable  ways.  In  the  derivation  of  the  expression  for  subclutter  visibility,  it  will  be  as¬ 
sumed  that  these  problems  can  be  neglected. 

The  dynamic  range  of  the  entire  system  is  then  limited  by  the  A/D  converters  and  the  process¬ 
ing  gain  following  them.  The  attenuation  control  must  be  adjusted  as  a  function  of  range,  for 
each  azimuth,  so  that  the  peak  signal  levels  are  just  below  the  A/D's  maximum  input  level.  In 
a  system  where  the  clutter  return  is  much  greater  than  the  target  and  receiver  noise  signals,  the 
clutter  peaks  should  not  be  allowed  to  saturate  the  A/Ds.  If  they  do,  the  clutter  spectrum  will 
broaden  out,  the  time  sidelobes  of  the  compressed  pulse  will  deteriorate,  and  a  loss  in  detection 
performance  will  result. 

Since  phase  coding  is  used,  the  attenuator  must  be  held  constant  over  the  basic  subpulse 
duration.  As  gain  changes  are  made  across  the  received  signal,  corresponding  gain  corrections 
must  be  introduced  in  the  digital  transversal  filter  so  that  the  pulse  compression  is  correctly 
performed.  The  simplest  gain  change  that  can  be  made  digitally  is  a  6-dB  change,  which  corre¬ 
sponds  to  a  word  shift.  Based  on  the  ground  clutter  seen  with  the  LRDR,  the  front-end  receiver 


62 


gain  should  not  be  monotonically  increasing  with  range.  One  should  be  ahle  to  accommodate  clut¬ 
ter  hot  spots  without  a  loss  of  subclutter  visibility  in  range  cells  following  the  range  cell  having 
large  clutter. 

The  second  problem  that  arises  with  the  use  of  an  A/D  is  the  nonlinearity  caused  by  the  quan¬ 
tization.  Signals  whose  peak-to-peak  excursions  are  less  than  a  quantizing  level  can  be  badly 
distorted  or  even  lost  unless  an  effective  dither  device  is  used.  This  can  he  accomplished  by 
injecting  receiver  noise  at  the  A/D  input  with  an  rms  value  of  at  least  one-half  a  quantizing  level 
which  is  independent  from  sample-to-sample.  This  is  done  by  making  the  second  handpass-filter 
bandwidth  equal  to  the  sampling  rate  of  the  A/D  and  by  adjusting  the  gain  of  the  second  RF 
amplifier. 

The  bandwidth  of  the  second  RF  amplifier  must  he  greater  than  or  equal  to  the  bandwidth  of 
the  basic  subpulse  transmitted.  If  it  is  equal  to  it,  then  one  has  essentially  a  filter  matched  to 
the  basic  subpulse,  and  the  digital  transversal  filter  consists  of  (N  -4)  shift  registers  and  adders. 
However,  the  bandwidth  of  the  second  RF  amplifier’  may  be  greater  than  that  of  the  basic  sub¬ 
pulse,  in  which  case  the  filter  matched  to  it  will  be  a  hybrid  filter.  If  the  bandwidth  of  the  filter 
is  twice  that  of  the  subpulse  then  the  transversal  filter  will  consist  of  2N  -  1  shift  registers  and 
adders . 

It  will  be  shown  that  the  subclutter  visibility  can  be  increased  by  opening  up  the  bandwidth 
of  the  second  bandpass  filter,  adjusting  the  video  noise  level,  and  increasing  the  A/D  sampling 
rate.  However,  the  receiver  eventually  will  become  susceptihle  to  jamming  and  inter  ference 
if  the  bandwidth  hecomcs  too  great. 

In  order  to  derive  the  optimum  receiver  structure,  a  complete  statistical  descr  iption  of  the 
clutter,  noise,  and  target  signals  is  necessary.  In  most  cases  of  interest,  the  received  clutter 
power  is  much  greater  than  the  target  power.  Hence,  unless  their  spectra  arc  separated  in  fre¬ 
quency,  reliable  detection  is  impossible  due  to  the  limited  time-on-target. 

It  will  be  assumed  that  the  clutter  signal  from  one  range  cell  can  be  modeled  as  a  zero  mean 
Gaussian  random  pr  ocess  whose  statistics  change  slowly  compared  with  the  processing  time. 

The  process  is  completely  characterized  by  its  power-  spectral  density  which  is  known  except  for 
a  finite  set  of  parameters.  An  example  would  be  a  spectra  of  known  shape  but  unknown  gain. 

The  receiver  noise  can  be  modeled  as  a  zero  mean,  band-limited,  white  Gaussian  random 
process  with  at  most  an  unknown  spectral  height.  This  spectral  height  will  he  a  function  of  the 
gain  control  which,  in  turn,  will  depend  on  the  clutter  power.  However,  for  a  well-designed 
system,  the  noise  spectral  height  will  remain  constant  over-  most  of  the  operating  range  of  the 
gain  control. 

The  target  returns  may  be  either  deterministic  or  random.  For  aircraft  and  ground  targets 
in  open  areas,  where  their  velocity  is  constant  over  the  observation  time,  the  target  signal  can 
be  modeled  as  a  known  signal  with  several  unknown  parameters.  Generally,  these  would  be  the 
amplitude,  phase,  arrival  time,  and  doppler  of  the  received  signal.  If  a  rotating  antenna  is  used, 
the  known  form  of  the  signal  will  depend  on  the  two-way  antenna  pattern.  If  a  phased  array  is 
used  in  a  stepped  scan  mode,  the  known  signal  form  will  depend  on  the  dwell  time  and  the  com¬ 
pressed  pulse  shape. 

The  return  of  a  target  in  foliage  will  be  random.  The  multipath  within  the  forest  causes 
the  field  to  vary  spatially  across  the  range  cell  which,  in  turn,  modulates  the  target's  return 
as  it  passes  through  the  cell.  For  a  constant-velocity  target,  its  spectrum  may  be  expressed 
as  a  known  function  with  a  finite  set  of  unknown  parameters.  These  parameters  could  be  a  gain. 


63 


a  bandwidth  measure,  and  a  frequency  shift.  It  is  not  clear  what  the  statistics  of  the  target  signal 
are.  However,  if  a  great  deal  of  multipath  were  present,  it  would  seem  reasonable  to  model  it 
as  a  Gaussian  random  process  with  unknown  parameters.  This  process  must  be  time-gated  by 
a  pulse  which  will  depend  on  the  compressed  pulse  and  the  target's  velocity. 

In  the  case  of  a  deterministic  signal,  the  receiver  consists  of  a  whitening  filter  followed  by 
a  bank  of  filters  matched  to  frequency-shifted  versions  of  the  target  signal  that  appears  at  the 
output  of  the  whitening  filter.  The  matched  filters  are  followed  by  square  law  envelope  detectors. 
For  an  arbitrary  observation  time,  the  whitening  filter  is  usually  time  varying.  However,  if  the 
observation  time  is  several  times  the  clutter  correlation  time,  the  filter  becomes  time  invariant. 
The  output  of  the  whitening  filter  is  white  noise  when  its  input  is  clutter  plus  receiver  noise. 

In  the  case  of  a  random  signal  which  can  be  modeled  as  a  Gaussian  process,  the  receiver 
can  be  implemented  in  the  following  manner  provided  that  the  observation  time  is  several  times 
the  clutter  and  target  correlation  time.  The  receiver  consists  of  a  whitening  filter  followed  by 
a  bank  of  target  filters.  The  outputs  of  the  target  filters  are  square  law  envelope  detected  and 
integrated  for  the  corresponding  time-on-target.  In  general,  the  target  filter  will  depend  on  the 
clutter,  receiver  noise,  and  target  spectra.  However,  the  target  filter  has  a  simple  form  if 
over  the  bandwidth  of  the  target  spectrum  the  target  power  is  appreciably  greater  than  the  clutter- 
plus -noise  power.  The  target  filter  can  then  be  approximated  as  having  unity  gain  over  this  band¬ 
width,  and  zero  gain  elsewhere. 

Unless  the  time -bandwidth  product  of  the  target  signal  is  large,  in  order  to  obtain  reliable 
detection,  the  target  power  must  be  larger  than  the  cluttcr-plus-receiver-noise  power  over  the 
target's  bandwidth. 

Given  a  probability  of  false  alarm,  the  detection  probability  is  directly  a  function  of  the  re¬ 
ceiver's  output  signal-to-clutter-plus-noise  ratio  when  the  signal  is  deterministic.  Unfortunately, 
for  the  random  signal  case,  the  output  signal -to-noise  (S/N)  ratio  cannot  be  related  directly  to 
the  detection  performance.  However,  by  adding  an  incoherent  integration  loss  factor,  the  signal- 
to-elutter-plus -noise  ratio  may  be  used  as  a  criterion  for  the  detection  performance.  This  loss 
factor  will  be  a  function  of  the  target  signal's  time-bandwridth  product,  the  probability  of  false 
alarm,  and  the  probability  of  detection. 

In  most  cases  of  interest,  the  target  signal  will  be  competing  only  with  receiver  noise  once 
proper  filtering  has  taken  place.  The  clutter  return  has  set  the  receiver  gain  level  so  that  the 
A/Ds  are  not  overloaded.  The  receiver  noise  level  is  set  by  the  linearity  requirement  of  the 
A/D.  Under  these  assumptions  and  by  using  the  optimum  receiver,  the  output  S/N  ratio  can  be 
expressed  as 


(S)  . 

lN  o 


in  which 


s 

(^)i  =  received  peak  signal  power/received  clutter 
1  power  from  compressed  range  cell 

G  =  processing  gain 

L  =  losses 

The  processing  gain  is  expressed  as  follows 
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in  which 


G  =  N  •  T  •  prf  •  2 
ot 


•  B  J  B 

ad'  p 


2(b-l ) 


N  =  length  of  phase  code 
TQt  =  time  on  target 


prf  =  pulse  repetition  frequency 
b  =  number  of  bits  in  A/D 
B  ^  =  A/D  sampling  rate 


-  bandwidth  of  compressed  pulse. 
The  loss  factor  is  expressed  as 


"ad 


L. 

1 


The  faetor  represents  the  loss  due  to  range  extended  clutter  when  pulse  compression  is 
used.  The  A/D  must  accommodate  the  uncompressed  clutter  return  without  limiting.  It  follows 
that 


Lc  =  i  for  point  clutter 

-  N  for  range  extended  clutter  with  a  uniform  cross  section 
over  the  code  duration. 


For  the  most  part,  the  loss  Lc  will  be  between  1  and  N. 

The  factor  represents  the  loss  due  to  an  increase  in  the  noise  level  at  the  A/D  as  the 
receiver  gain  is  increased  and  the  front-end  noise  predominates.  For  a  well-designed  system, 
this  will  only  oeeur  at  the  maximum  range. 

The  factor  represents  the  loss  due  to  the  A/D  limiting  and  linearity  margins  and  is  given 


,d  ■ 


in  which 


6  =  A/D  limiting  margin.  The  rms  clutter  voltage  is  set  to  E/<5,  where 
E  is  the  magnitude  of  the  A/D's  maximum  input. 

y  -  A/D  linearity  margin.  The  rms  noise  voltage  is  set  to  Q/y,  where  Q 
is  the  A/D‘s  quantizing  increment. 

The  factor  L.  represents  the  incoherent  integration  loss  for  random  target  signals.  If  the 
target  signal  is  deterministic,  then  L.  =  i;  otherwise. 


in  which 


(T 


ot 


V 


T  t  =  time  on  target 
Bg  =  target  signal  bandwidth  >  t/T^. 


65 


The  factor  a  is  such  that  0  a  l/2  and  is  a  function  of  the  coherent  S/N  ratio,  probabilities 

of  false  alarm  and  detection,  and  the  time-bandwidth  product  T  ,  •  B  . 

ot  s 

In  a  forest  where  the  field  can  be  considered  a  spatially  homogeneous  random  process  with 
a  correlation  distance  of  r*c,  it  follows  that  the  incoherent  integration  loss  can  be  expressed  as 


in  which 

r^  -  correlation  distance  of  the  field  across  the  range  cell 

AH  effective  range-eell  length  =  (C/2)  where  T  n 

is  the  effective  compressed  pulse  duration. 


In  this  case,  the  target  spectrum  possesses  a  constant  Q  and  as  the  target's  velocity  increases, 
its  bandwidth  increases  while  the  Tof  decreases.  It  has  been  assumed  that  the  range  cell  defines 
the  Tot  and  that  the  target  has  a  constant  velocity. 

Given  a  desired  output  S/N  ratio  (S/N)  the  subclutter  visibility  (S/C)^  of  a  radar  with  a 
processing  gain  G  and  a  loss  L  can  be  defined  as 

KC9v  '  N'od  G  * 


Then,  so  long  as  (S/C).  MS/C)  ,  the  radar  will  perform  adequately. 

Given  that  the  clutter  area  cross  section  is  uniform  across  the  range-azimuth  cell,  the  input 
signal-to-clutter-power  ratio  can  be  expressed  as 


(|»f  -  - 


n  '  Ae  *  ar 


in  which 

cro  =  clutter-area  cross  section 
Fc  -  tw*o-wray  clutter  propagation  factor 
<rt  -  target  cross  section 
Ft  =  two-way  target  propagation  factor 
H  range  to  range-azimuth  cell 
AO  two-way  antenna  beamwidth 
AR  -  effective  range-eell  length  =  (C/2)  T  ^ 
Tc  ii  =  effective  compressed  pulse  duration. 


If  all  the  parameters  above  are  known,  it  then  can  be  determined  whether  reliable  detections  are 
possible. 

In  many  cases,  one  would  like  to  compare  different  radars  to  determine  which  performs  best. 
Thus,  the  comparison  should  be  made  on  an  output  S/N  ratio  basis.  Assume  that  range  extended 
clutter  predominates  so  that  -  N,  that  the  A/D  noise  level  is  such  that  =  1,  and  that  the 
A/D  limiting  and  linearity  margins  are  set  so  that  L  ^  =  1.  Then  the  output  S/N  power  ratio  can 
be  expressed  as 

,s.  gt  ’  *  t  _  _ Tot  ’  prf _  .  R  .  _2(b- 1 ) 

N  °  ao  ‘  Fc  B  •  (T  B  f  •  R  •  Ae  •  AR  ad 
p  ot  s 
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Let  us  now  consider  two  cases. 


Case  \ 

Let  the  dwell  time  on  a  range-azimuth  cell  be  long  so  that  the  TQt  depends  on  the  range-cell 
length  and  target  velocity.  Then,  T  =  AR/V  where  V  is  the  magnitude  of  the  target's  radial 
velocity.  For  the  case  of  a  coherent  target  signal,  it  follows  that 

,S.  _  gt  •  Ft  prf  '  Tcell  2(b-1 ) 

Vo  <T  •  F  '  R  •  AO  •  V  ad  '  z 
o  c 

It  should  be  remembered  that  B  _  >.  l/T 

ad  '  cell 

In  comparing  different  radar  systems,  it  follows  that  (at  *  F^)/(aQ  •  F  )  is  a  medium- 
dependent  parameter  while  the  remaining  terms  are  system  dependent. 

Provided  that  the  target's  velocity  is  constant  over  the  range  cell  and  that  the  A/D  sampling 
rate  remains  constant,  the  output  S/N  ratio  is  proportional  to  the  range-cell  length.  There  arc 
system  constraints  such  as  noise  level  and  gain  adjustments,  the  RF  bandwidth  of  the  receiver, 
and  the  complexity  of  the  digital  transversal  filter  that  will  limit  the  usefulness  of  choosing 
Bad  >  VTcen. 

If  the  A/D  handwidth  is  matched  to  the  compressed  pulse  bandwidth  so  that  T  ^  1, 

it  then  follows  that  the  output  S/N  power  ratio  docs  not  depend  on  the  range -cell  duration 

There  is  a  tradeoff  between  the  A/D  sampling  rate  and  the  number  of  hits,  provided  that  the 
range-cell  duration  can  be  chosen  so  that  T  ^  ^  l/B^.  In  this  case,  a  useful  figurc-of-merit 
of  the  A/D  can  be  defined  as 

M  101o*iofBad-  22<b‘ni  • 

In  this  manner,  different  A/Ds  may  be  compared  when  designing  a  digital  MT]  radar. 

In  the  case  of  a  target  in  foliage  where  the  target  return  is  random,  the  coherent  S/N  power 
ratio  must  be  multiplied  by  the  factor  (AR/rc)  a  -  (CTc  n/2r  )  a .  Hence,  the  output  S/N  ratio 
will  be  proportional  to  T  jj,  Sinee  0^  a  ^  1/2,  the  system  will  perform  better  as  T  ^  is 
increased,  provided  that  Ba^  ^  l/T^e 


Case  2: 

Let  the  dwell  time  on  the  range-azimuth  cell  be  controlled  by  the  antenna  scan  time  and  not 
the  target  velocity.  Then,  the  output  S/N  ratio  for  a  coherent  target  has  the  following  form; 

,S,  fft  '  ?t  2  •  Prf  •  Tot  #  2(b-l) 

Vo  “  a  •  F  ’  R  •  AG  •  C  ad 
o  e 


Again  it  is  seen  that  this  expression  can  be  separated  into  medium  and  system  parameters. 

For  a  fixed  A/D  sampling  rate  and  provided  that  ^  **  is  seen  that  the  output 

S/N  power  ratio  does  not  depend  on  the  range-cell  duration.  This  property  is  attained  at  the 
expense  of  the  increased  digital  hardware  in  the  transversal  filter.  If  the  A/D  sampling  rate  is 
matched  to  the  compressed  pulse  bandwidth,  then  the  output  S/N  ratio  is  inversely  proportional 
to  the  range-cell  duration  T  jj, 

As  in  Case  1,  there  is  a  tradeoff  between  the  A/D  sampling  rate  and  number  of  bits.  The 
criterion  factor  M  ean  be  used  to  compare  different  A/Ds. 
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CLOCK  SIGNAL  TO  ALL 
OIGITAL  LOGIC 


Fig.  V-l.  Analog  section  of  receiver. 


Fig.  V-2.  Digital  section  of  receiver. 
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